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ABSTRACT 


This thesis provides a crosstalk analysis of optical chip interconnects via single- 





mode waveguides with synchronous transmission and asynchronous transmission. This 
crosstalk model is genera] and can be used for any type of waveguide network. Three 
cases of laser sources will be considered: (1) each channel operates with an independent 
laser sources, (2) all laser sources have the same mean wavelength but with different 
phase noise processes, and (3) all laser sources are identical with the exception of the 
initial phases. The analysis takes into account the coupling-induced crosstalks between 
adjacent waveguides, the laser linewidth, the shot noise, the dark current generated by 
the photodiode, and the post-detection thermal noise. Bit error probabilities versus 


received peak powers are presented together with power penalties. 
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I. INTRODUCTION 


This thesis considers the problem of crosstalk impact in an optical interconnect 
System using a single-mode waveguide network [Ref. 1-Ref. 8]. An optical interconnect 
consists of two or more terminal nodes which are interconnected hv a single optical path 
or a high-density parallel network which preserves the parallel nature of the data 
generated at the nodes [Ref. 9]. Such interconnect systems would consist of chip-to-chip 
or board-to-board interconnections [Ref. 8]. The use of waveguides provides a potential 
of integrating the entire system of transmitter and receiver on the same substrate [Ref. 


10). 


Figure 1: Schematic diagram of an optical interconnection transmission system 
architecture [ Ref. 9). 


A. BACKGROUND 

The success of optical communication has accelerated research on high capacity data 
handling systems. It is expected that the monolithic integration of optical and electronic 
components on the same chip will ultimately lead to ultrahigh-speed, high sensitivity, 
compactness, reliability, low cost, as well as passive and active integrated optic 
components [Ref. 11]. 

Integrated optics refers to the implementation of various functions with light such as 
modulation, switching, generation, and detection in an optical guided wave structure 


formed on a substrate. Use of the term "integrated" implies implementation of many of 








these functions on the same substrate. Integrated optical devices are characterized by 


many advantages associated with lightwave technology in general, namely, larger 
information capacities than with electronic technology are possible, electromagnetic 
interference is not present, and parallel processing is possible. Many of the advantages of 
electronic integrated circuits, including the potential for fabrication economy and 
reliability of devices combined onto one substrate apply to integrated optics. Vibration 
problems associated with bulk optical experiments are eliminated when devices can be 
integrated onto one substrate space. [Ref. 12] 

Closely associated with integrated optics is the field of integrated optoelectronics 
which encompasses device structures referred to as optoelectronic integrated circuits 
(OEICs) [Ref. 12]. OEICs represent a device technology with potential to meet a broad 
range of future telecommunication and computing system needs. Optoelectronic 
integrated circuits are circuits that monolithically integrate optical and electrical 
components on a single semiconductor chip [Ref. 10]. Compound semiconductor 
materials are used to form optoelectronic integrated circuits because photonic devices 
such as semiconductor lasers, detectors, high-speed electro-optic modulators and 
switches, as well as quantum well waveguide devices, and high-speed electronic devices 
can be formed with the same material alloys [Ref. 12). 

High performance is now being achieved in devices using both GaAs and InP 
material systems [Ref. 9}. It is in the combination of photonic components and electronic 
circuitry that the OEIC gains a usefulness over and above that which can be obtained by 
placing non-monolithic circuits together in a package. For example, it has been asserted 
that only through monolithic integration can one fabricate extremely high-bandwidth 
transmitters or high sensitivity receivers [Ref. 9]. 

OEICs currently fall into three categories: interconnects, communications, and 
computing and signal processing. As mentioned before, an optical interconnect (Fig. 1) 


consists of two or more terminal nodes which are interconnected by a single, high-band 





optical path, or alternately, a high-density parallel optical network which might preserve 
the parallel nature of the data generated at the terminals. Such an architecture is useful 
for interconnecting two computer mainframes or, on a smaller scale, would consist of 
board-to-board or even chip-to-chip interconnections. The desirability of such a system 
is its potential for the rapid transfer of parallel data from one system to another without 
the concomitant problems of electromagnetic interference ari signal dispersion 
commonly observed in a high-bit-rate electrical interconnects. [Ref. 9] 

Fos example, as computing environments evolve toward the type of distributed 
network in Fig. 2 with data processing and database sharing among remote locations, 
Overall operation efficiency increasingly relies on the efficiercy of interconnection links. 
Optical interconnects are a new approach that could be used to achieve high-bandwidth 
low-loss interconnects for these applications. (Ref. 9] 

Figure 3 illustrates the functional charac.eristics of a GaAs computer interface chip . 
For insertion into a practical computer system, this circuit, for any OEIC component, has 
been designed to meet the following requirements (Ref. 10]: 

1) The circuits must be high speed, capable of multi-gigabits per second, with 
relatively high complexity (10K transistors/chip). 

2) They must be capable of operation in noisy envircaments; P,, < 10-15 in the 
presence of noise levels up to 100 mV. 

3) They must have high reliability, i.e., a failure rate <0.01% per kilohours for the 
entire link, for over 105 h at 50° C. 

4) They must have redundancy in critical paths with transmission error correction 
capacity. 

5) They must be compatible with existing computer interface technologies, 


packaging, power supply requirements, and silicon interface. 








Work 
Station 


Processor 





Figure 2: Distributed computer network [Ref. 10}. 


The advantages of OEICs are currently being exploited by the Department of the 
Defense. OEICs are being developed for such applications as neural networks, high 
speed signal processing, high speed communications, and advanced antenna systems. For 
example, the development of efficient microwave signal distribution by optical fiber 
creates the possibility of achieving unconventional antennas such as "smart skins". In 


this example, individual radiated elements would conform to the contour of an aircraft 





and the "antennas" directional properties would be determined by contro] of phase 
emission from different parts of the aircraft. Consequently, the DOD funds on-going 
programs developing high frequency components, such as high speed GaAs circuitry. 
This program is funded for seven years at a level of $500 million. (Ref. 10} 





Figure 3: OEIC Functions (Ref. 10] 


B. APPROACH 

This thesis considers a mathematical approach that includes nivie than two 
interfering channels. The modulation format used in this analysis is on-off keying (OOK) 
with direct detection. Two, three, and four channels were modeled using synchronous bit 
transmission. For comparison, two adjacent channels were modeled using asynchronous 
bit transmission. The emphasis of this thesis will be on the impact of coupling induced 
crosstalk between adjacent waveguides in a waveguide network. Furthermore, additional 


parameters that will affect the performance of the interconnect system are the laser 








linewidth, the shot noise and the dark current generated by the photodiode, and the post- 
detection thermal noise. To be general, the receiver will be considered to be an integrate- 
and-dump filter with integration time T where T is the bit time. Also, the current spectral 
density of the post-detection thermal noise is denoted as No. The spectral density No can 
be easily computed for a given low noise amplifier type given the effective noise 
temperature and the matched resistance load. 

The optical interconnect system is given in Fig. 4(a) for OOK direct detection. The 
waveguide network consists of many single-mode waveguides. It can be a planar array 
with uniform waveguide separation or any other structure (Fig. 4(b)). The envelope of 
the output lightwave of a given waveguide is detected by the photodiode, which also 
generates shot noise and dark current. The output of the photodiode is further corrupted 
by the amplifier thermal noise. The total signal plus crosstalk, shot noise, dark current, 
and thermal noise is integrated over one bit time T and the resulting bit energy at the end 
of each integration time is determined by the slicer to be either bit one or bit zero. It is 
assumed that the waveguide bandwidth is much larger than the bit rate and the signal 
spectrum. 

The resulting model for synchronous and asynchronous transmission was validated 
using MATLAB. The resulting data from the MATLAB programs was graphed as a 
series of performance and optimal threshold curves versus peak powers. These resulting 
curves were then compared with Ref. 8 and Ref. 13 to determine the credibility of the 
results. Appendices C through F contain sample programs for synchronous and 


asynchronous transmission. 

















Figure 4(a): Optical Interconnect System for OOK Direct Detection 


Planar Waveguides 


Figure 4(b): Planar Waveguides 


C. THESIS ORGANIZATION 

This thesis is organized into four chapters. Chapter I provides a background to 
optical chip interconnects. Chapter II provides the theoretical framework for the 
performance analysis and discusses the numerical results when the bit transmission is 
assumed to be synchronous where bit streams in all channels are time-aligned. Chapter 
III provides the theoretical framework for the performance analysis and discusses the 
numerical results when the bit transmission is assumed to be asynchronous where bit 
streams in all channels are not time aligned. As in Ref. 8, three cases of laser sources 


were considered and modeled. In case one, the system was modeled utilizing non- 











coinciding and uncorrelated light sources (each channel operates with an independent 
laser source). In case two, the system was modeled utilizing coinciding and uncorrelated 
channel] light sources (all laser sources have the same mean wavelength but with different 
phase noise processes). Finally, in case three, the system was modeled utilizing 
coinciding and correlated channel] light sources (all laser sources are identical with the 
exception of the initial phases). Chapter IV is a summary of the study. Appendices A 
and B are derivations of the conditional mean square value of X; for synchronous and 
asynchronous transmission. Appendices C through F contain sample MATLAB 
programs. Appendices C and D are sample programs for two adjacent channels with 
synchronous and asynchronous transmission. Appendices E and F are sample MATLAB 


programs for three and four adjacent channels with synchronous transmission. 








Il. SYNCHRONOUS OPTICAL CHIP INTERCONNECTS 


A. ANALYSIS 

For mathematical convenience, the complex envelope notation of a real signal is 
adopted. Thus, for a given transmitted bit b;9 of a given channel 0 whose laser phase 
noise process is Qo(t), the signal at the input of the photodiode is designated as 


; M 
s,(t) = Fe digeho + YPo gel rorrrn, O<tsT (1) 
k=] 


where A is the OOK signal amplitude, 6, (t) is the laser phase noise process, @, and $, 
are the frequency spacing and the initial phase difference between channels k and 0, 
respectively. Also, b;, represents the transmitted bit of channel k, and C, represents the 
coupling from channel k to channel 0. The summation term in equation (1) thus 
Tepresents the crosstalk from M adjacent channels into channel 0. In equation (1), all 
waveguides are assumed to have the same attenuation. For the case of M = 2, that is, 
when only the two nearest adjacent channels are considered, equation (1) reduces to the 
result in equation (2.4) in Ref. 8. In this case C, = B/ ¥1-2B? , k=1,2 where B is the 
total power coupled from a waveguide into its adjacent waveguide. In practice b,, k = 
0,1,---,M are not necessarily equal to zero or one. Let r,, k = 0,1,---,M be the extinction 


ratio of the laser of channel k defined as the ratio of the transmitted power of the logical 
zero to that of a logical one. Then by = a /(1+%) for logical zero and 


by, = Ji/fi+n) for logical one. Hereafter, the r,'s are assumed to be identical for all 
channels. 

Let R = n,q/hf be the photodiode responsitivity [Ref. 14 - Ref. 15] where n, Ss 1 is 
the quantum efficiency, q is the electron charge (1.6x10-!9C), h is Planck's constant 


(6.626 x 10-34 J-s), and f is the frequency. The output current of the photodiode is 


R\s;(t) +w,(t)+ Wg, (t) where w,(t) is the shot noise generated by the photodiode and 














Wa, (t) is the dark current noise. The output of the photodiode plus the post-detection 
thermal noise n(t) is integrated by the integrate-and dump filter with the normalization 


constant R resulting in the following decision variable Y; 


1 ¢T 2 1 pT 1 pT 1 fT 
Y; = = ls Ris, (1) dt+ Fo walt + he Wa (t)dt + al n(t)dt (2a) 
= X; + W, + Wa +N 
where 
2 
X= [sof at (2b) 
1 pT 
Wee j pWe(t)dt (2c) 
1 fT 
Wa =5 J, Wax (tat (2d) 
1G? 
N= J, n(tht. (2e) 


Since n(t) is a zero mean Gaussian process with spectral density No, the Gaussian random 


variable N also has zero mean and its variance oO}, is given by 


™ 
On = a (3) 


On the other hand, the shot noise w,\:) is a non-stationary process since the envelope of 
the signal at the input of the photodiode, namely |s;/t}, is time-dependent. Therefore, the 
shot noise w,(t) can be modeled as a zero mean wide-sense stationary Gaussian process 
whose spectral density Wo( bi), given a bit pattern bj = (b;o, bj1,-~-,bx), is proportional to 


the conditional mean of the squared envelope of the input signal. In other words, 
E{\s,(ey"[bs} approximates |s;(t)| over a bit time T. Based on this approximat'n, the 


10 














shot noise spectral density can be obtained as follows (Ref. 17] 
Wo(bi) = GRE |{s;(1))"[bi} 


1 M 

~ }ona{vf + Sted} Y 
k=) 

From equation (2c) and equation (4), the conditional variance of W,, given a bit pattern 


bi, is 


— \ TWo(dbi 
of, (bi) = Pe, (5) 


The dark current noise spectral density function is qla, where Igy is the dark current. The 
variance of the dark current noise is Ow, = Tqly / R?. The random variable Xj in 
equation (2b) consists of the signal] term, the signal-crosstalk term, and the crosstalk- 


crosstalk term. Substituting equation (1) into equation (2b), the random variable becomes 
1 a2 1S, 42 ~I0, {7 {09 (1)-0, (1-0, t] 
Xi=5A Thio+5 LCA biodiy {e7 J eAol-OO-ovtay 
k=1 
+eite in eo Seotreora]} (6) 


+S > sC,C Abb, eX%e-%) fp el (1)-O¢(1)+(@,~¢)t] gy 
k=1f=1 
The last two terms of the expression of Xj; represent the crosstalk in channel 0. The 
statistics of these terms are extremely difficult to obtain (if possible). Consequently, for 
tractable analysis, Yj; is modeled as a Gaussian random variable. Such a Gaussian 
approximation has also been used in Ref. 8 with a different mathematical approach. 
Gaussian approximations are commonly used to obtain the bit error probability for 
lightwave systems when the exact statistics of the decision variable cannot be analytically 
obtained [Refs. 13, 17 -18]. From equation (6), the mean of X; conditional on a given bit 


pattern bj is given by 
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X , (b;) = E{X;[b;} = SA Tin +5 LCA TOR. (7) 
k=1 


From equation (7) and equation (A5) of Appendix A, the conditional variance OX, (bi) of 
X; can be calculated as follows: 
M 
ox (bi) = LCLA‘T?bjbk {2xv/4n2(v? +82) 
+82 


aa ave?" (225, sin 2nd, — 2RVCOS 25, ) (8) 


+2n5,¢-2™ (205, cos2n5, + 2xvsin2n8,)+ 4n°v? ~ 4n°5i]} 


where v = BT and 5, = @,T/22. The parameter B is the laser linewidth. In summary, 
the Gaussian approximation allows the decision variable Y; in equation (2) to be 
considered as a Gaussian random variable with conditional mean ¥i(bi) = Xi(b i) in 


equation (7) and conditional variance o¥, (Bi) given by 


of, (5) = of (5) +04, (5)+ 0% toh. ) 


For a threshold «, the conditional bit error probability given bit patterns bP = (bo9, 
bi bim) and b = (dio, Disses bim) is [Ref. 19] 


P, (bit»---»bim) = =Po(br) +58 (6%) (10a) 





¥,(bi)-a 


~1 us 1 
(Bi) 2 erfe Voy, (6) (e) (10c) 
and erfc( - ) is defined as 
erfc(a) = fre ex. (10d) 


The bit error probability P,, is obtained by taking the expectation of P,,(b,;,---*b4) with 
respect to the bit patterns (b;,,---, by). Since there are 2M such patterns, the bit error 


probability becomes 


1 
P, = pwd (biyeeabag) Pb(Dit>---»Dim ) (11) 


where the summation is over all 2M patterns (b,;,bj,2,--.6y4)- The optimal threshold that 


minimizes the bit error probability is the value that satisfies equation (12). 





(12) 





Equation (12) is obtained by setting OP, /da to zero. 
In the case when all laser sources have the same mean wavelength but are 
uncorrelated, the above results apply by setting 5, = 0 with k = 1, 2,---, M in equation 


(8). Furthermore, when all laser sources are identical such that all sources have the same 
wavelength and phases noise process except for the random initial phases, then OX, (bi) 
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in equation (8) reduces to 


1 M 
oX, (bi) = 5 LAT’ babi (13) 
k=] 


B. NUMERICAL ANALYSIS 

In this section numerical results are presented for a system with a bit rate of 500 
Mb/s. The responsitivity of the photodiode is taken to be 0.5 and the laser extinction 
ratio is 1/20. The dark current I4,=10 nA. The effective noise temperature for a low 
noise amplifier-integrate and dump-slicer receiver is 180 K. Furthermore assuming a 
matched load of R, = 50 Q, the post detection thermal noise current spectral density is 
No = 2kTo/R, = 10°22 A2/Hz [Ref. 15 and Ref. 20) where k = 1.38 x 10°23 J/K is the 
Boltzmann's constant. 

This section is broken down into the three cases of laser sources and the power 
penalty plot summaries. The data collected are the results of the MATLAB model 
(Appendix C) of the optical chip interconnect for synchronous transmission. The first 
case of laser source occurs when all channels operate with independent laser sources. 
Case two occurs when all laser sources have the same mean wavelength but have 
different noise processes. Finally, in case three, all laser sources are identical with the 
exception of the initial phase. These three cases are summarized in the three power 


penalty plot summaries for a bit error probability Pp of 10-15. 


1. Casel 
Figures 5-6 show the bit error probability P, versus the received peak power 
A2/2 for various levels crosstalk from two adjacent channels relative to that of a single 


channel operation (zero crosstalk). The laser linewidth-bit rate ratio and frequency 
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spacing-bit rate ratios (v,5,,5,) were taken to be (0.1,0.3,0.3) and (0.1,0.7,0.7), 
respectively. As seen in Fig. 5, the bit error rate floor exists around 10-!2 for —-20 dB 


crosstalk irrespective of the received peak power. Crosstalk levels must be less than —26 
dB for a power penalty of 1 dB or less at Py = 10°}5. 


: 
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Peak Power (dBW) 





Figure 5: Py vs. A?/2 as a function of crosstalk levels with (v,5,,5) = (0.1,0.3,0.3). 


From Fig. 6, it is seen that by increasing the frequency spacing, the crosstalk 
level can be reduced to less than —23 dB for 1 dB or less in power penalty at P, = 10°18. 
Figures 7-8 show the normalized optimum threshold versus the received peak power as a 


function of crosstalk levels. The optimum threshold decreases with increasing crosstalk. 
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Figure 7: Normalized optimal threshold « vs. A2/2 as a function of crosstalk levels 
with (v, 54 52) = (0.1,0.3,0.3). 
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Figure 8: Normalized optimal threshold o vs. A2/2 as a function of crosstalk levels 
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Figure 9: 


with (v,51,59) = (0.1,0.7,0.7). 


Peak Power (dBW) 


P) vs. A?/2 as a function of crosstalk levels with (v,51,5,) = (1,0.3,0.3). 
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Figure 10: P, vs. A2/2 as a function of crosstalk levels with (v,5,,5) = (5,0.3,0.3). 


Figure 9 shows the results for (v,8,,52) = (1,0.3,0.3). It is seen that a laser 
with a larger laser linewidth-bit rate ratio improves the performance. For 1 dB or less in 
power penalty at P, = 10-15, the permitted crosstalk level is less than—23 dB instead of 
— 26 dB as in Fig. 5. When the normalized linewidth is increased tc 5 as in Fig. 10, there 
is less than 1 dB power penalty at crosstalk levels less than —17 dB. This happens 
because only a portion of the crosstalk energy fall within the detection bandwidth of 500 
MHz. This result encourages the use of lasers with a large linewidth as long as the 
waveguide bandwidth is larger than the signal spectrum. When the signal spectrum 
broadened by the laser phase noise approaches the waveguide bandwidth, loss in signal 
power begins to occur and performance deteriorates rapidly. 

For comparison, this model was expanded to include the cases for three and 


four adjacent channels for synchronous transmission. Figures 11 and 12 show the 








performance for the cases of three and four adjacent channels in a nonplanar waveguide 
network with (v,5,,52,53)=(0.1,0.3,0.3,0.3) and (v,8,,52,53,8,)=(0.1,0.3,0.3,0.3,0.3), 
respectively. For three adjacent channels, a crosstalk level of less than —29 dB is 
required for a power penalty of 1 dB or less at P, = 10°-!5, as compared to a crosstalk level 
of —26 dB for the case of two adjacent channels as shown in Fig. 5. The bit error floor 
for three adjacent channels exists at 10-9 for a ~20 dB crosstalk level irrespective of the 
teceived peak power. In Fig. 12 for four adjacent channels, a crosstalk level of less than 
approximately —30 dB is required for a power penalty of 1 dB or less at Pp = 10°15. 
Here, the bit error floor exist around 10-7 for —20 dB crosstalk level. The differences 
between Figs. 11 and 12 are not as dramatic as those seen between Figs. 11 and 5. 
Figures 13 and 14 show subtle changes in the normalized optimal threshold in 
comparison to Fig. 5. There appears to be a very slight change in the normalized optimal 


threshold as the number of adjacent channels increases. 
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Figure 11: P, vs. A2/2 as a function of crosstalk levels with (v,5,,5,,53) = 
(0.1,0.3,0.3,0.3) 


19 











~54 -52 -50 





Peak Power(dBW) 


Figure 12: P, vs. A2/2 as a function of crosstalk levels with (v,5,,52,53,54) = 
(0.1,0.3,0.3,0.3,0.3). 





Figure 13: Normalized optimal threshold o vs. A2/2 as a function of crosstalk levels 
with (v,5,,5,53) = (0.1,0.3,0.3,0.3) 
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Figure 14: Normalized optimal threshold a vs. A2/2 as a function of crosstalk levels 
with (v,5,,52,53,54) = (0.1,0.3,0.3,0.3,0.3). 


2. Case Il 
Figure 15 shows the performance of case two where all laser sources have the 
same mean wavelength but have different noise processes for (v,5,,52) = (0.1,0,0). A bit 
error rate floor exists at Py, = 10-19 for—20 dB crosstalk level. The permitted crosstalk 
level is less than —29 dB for a power penalty of 1 dB or less at P, = 10-15. In general the 
performance of case two is always worse than that of case one given the same v and 
crosstalk level. It is obvious that the crosstalk effect is reduced by using laser sources 


with different wavelengths. 


3. Case Ill 
Figure 16 shows the performance of case three when all laser sources are 
identical with the exception of the initial phases. The performance is slightly worse than 


that in Figure 15. 
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Figure 15: P, vs. A2/2 as a function of crosstalk levels with (v,53,52) = (0.1,0,0). 
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Figure 16: P, vs. A2/2 as a function of crosstalk levels with all laser sources having 
identical phase noise except for initial phases. 
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4. Power Penalty Plot Summary 

> results for two adjacent channels with synchronous transmission are 
summarized in Figs. 17, 18 and 19 as power penalty relative to a single-channel operation 
versus crosstalk levels, versus the frequency spacing-bit rate ratios 5, = 5, and versus 
the normalized linewidths, respectively. All results were taken at Pp = 1075. Figure 17 
shows that as the crosstalk level increases, so does the power penalty. It also shows that 
the smaller the normalized linewidth, the more sensitive the system is to smaller crosstalk 
levels. Figure 18 shows that as the power penalty increases, the normalized frequency 
spacing decreases. Finally, Fig. 19 supports the results in Fig. 17. For normalized 


linewidths less than 5, the power penalty dramatically increases. 


w lb requency, Spacing = 0.3. enebef ee cea mnmamnedimenanaf coset 





Figure 17: Power penalty vs. crosstalk level for a normalized frequency spacing of 
51 = 5, = 0.3. 
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Figure 18: Power penalty vs. normalized frequency spacing (5, = 52) for a 
normalized linewidth of v = 0.1. 





Figure 19: Power Penalty vs. normalized linewidth for 5; = 6, = 0.3. 











C. DISCUSSION 

This chapter presented a mathematical framework to analyze the performance of 
synchronous optical chip interconnects in terms of the bit error probability versus the 
received power as a function of the crosstalk level, frequency spacing, and laser 
linewidth. This analysis can handle any number of adjacent channels. The conclusion 
drawn from this investigation is that adjacent channels must use laser sources of different 
wavelengths to reduce the effect of crosstalk. Laser sources with a large linewidth also 
help, as long as the waveguide bandwidth is much larger than the resulting signal 
spectrum. In fact, this is the only way to reduce the effect of a given crosstalk level when 
all laser sources have the same mean wavelength. When all laser sources are locked to 
the master source, the performance depends explicitly on a crosstalk level given a 
received peak power. 

Similar conclusions appear in Ref. 8 via a different mathematical approach that 
applies to two adjacent channels only. The mathematical framework presented in this 


chapter will be applied to the asynchronous case in the next chapter. 
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Ill. ASYNCHRONOUS OPTICAL CHIP INTERCONNECTS 


A. ANALYSIS 
As in the synchronous transmission analysis, for mathematical convenience, the 
complex envelope notation of a real signal is adopted. Thus, for a given transmitted bit 
b;o of a given channel 0 whose laser phase noise process is Q(t), the signal at the input 

of the photodiode is designated as 
s;(t) = * 


M 
yea el®s (+, (0404.1 ]fu(t)~ u(t - t,)] (14) 
a V2 


b; eo u(t) — u(t—T)] 


MCA, je 
+> ts bi ot ea, (+O.0]f y(t —%)-—u(t—- T)] 
k=1 


where A is the OOK signal amplitude; bj.) and bi, are the previous and the present bit 


in channel k relative to bit bio; 6, (t) is the laser phase noise process; and ®,, >, -; and 
x,0 are the frequency spacing and the initial phase differences between channels k and 
0, respectively. The parameter T, represents the uniformly distributed random delay 
between bio and b;x,9, The function u(t) denotes the unit step function. The parameter Cy 
represents the coupling from channel k to channel 0. The summation term in equation (1) 
thus represents the crosstalk from M adjacent channels into channel 0. In equation (1), all 
waveguides have the same attenuation. For the case of M = 2 and synchronous bit 
transmission, that is, only the two nearest adjacent channels are considered, equation (1) 
reduces to the result in equation (2.4) in Ref. 8. In this case C, = B/ 1~2B?, k=1,2 
where B is the total power coupled from a waveguide into its adjacent waveguide. In 
practice, bi, , and bi 9 k = 0,1,---,.M are not necessarily equal to zero or one. Let %, k = 
0,1,--:,M be the extinction ratio of the laser of channel k (defined as the ratio of the 


transmitted power of the logical zero to that of a logical one). Then 
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Dox.o»Dox.-1 = vi /(i+%,) for logical zero and by. 9,64. -1 = vi /(i+x,) for logical 
one. Hereafter, r,'s are assumed identical for all channels. 

Let R=n,q/hf be the photodiode responsitivity [Ref. 14 - Ref. 15] where n, $ 1 is the 
quantum efficiency, q is the electron charge (1.6x10-!9C), h is Planck's constant (6.626 x 
1034 J-s), and f is the frequency. The output of the photodiode is 
Ris,(t) + w,(t)+ Wa, (t) where w,(t) is the shot noise generated by the photodiode and 
Wax(t) is the dark current noise. The output of the photodiode plus the post-detection 
thermal noise n(t) is integrated by the integrate-and dump filter with the normalization 


constant R resulting in the decision variable Y; as follows: 


2 
Y= = fF Risi(ey dt+ = [owe + = i Wa, (t)dt + - JF nceyae 


(15a) 
= X; + W,; + Wa +N 
where 
T 2 

Xi = [i Is(e){ at (15b) 
1 (pT 

W, = ah w,(t)dt (15¢c) 
1, 

Wa, = Flo Wax (tat (15d) 
1 fT 

N= ai n(t)at. (15e) 


Since n(t) is a zero mean Gaussian process with spectral density No, the Gaussian random 


variable N also has zero mean and its variance on is given by 


of = = (16) 
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On the other hand, the shot noise w,(t) is a non-stationary process since the envelope of 
the signal at the input of the photodiode, namely Is;(t)l, is time-dependent. Therefore, the 
shot noise w,(t) can be modeled as a zero mean wide-sense stationary Gaussian process 
whose spectral density Wo( bi), given a bit pattern bi = (b;0,b)1,-15--"sDiM,-1,Di1,00"""DiM.0), is 
proportional to the conditional mean of the squared envelope of the input signal. In other 
words, E{[s:(t)? bs} approximates Is;(t)I2 over a bit time T. Based on this approximation, 
the shot noise spectral density can be obtained as follows [Ref. 17] 
om 2 — 
Wo(bi) = aRE{Ks,(1)"[6i} 


1a 1 
ara os > Cib5. 1+ 3 Sched} (17) 
k=] k=] 


at 
2 
From equation (15c) and equation (17), the conditional variance of W, given a bit pattern 
b; is 


— \ TWoldi 
0%, (6) = wo). (18) 


The dark current noise spectral density function is qla, where Ig, is the dark current. The 
variance of the dark current noise is OWa = Tqlay / R*. The random variable Xj in 
equation (15b) consists of the signal term, the signal-crosstalk terms, and the crosstalk- 
crosstalk terms. 

The statistics of the crosstalk terms are extremely difficult to obtain (if possible). 
Consequently, for tractable analysis, Y; is modeled as a Gaussian random variable. Such 
a Gaussian approximation has also been used in Ref. 8 for synchronous transmission with 
a different mathematical approach. Gaussian approximations are commonly used to 
obtain the bit error probability for lightwave systems when the exact statistics of the 
decision variable cannot be analytically obtained [Refs. 13,17 -18]. From equation (15b), 


the mean of Xj conditional on a given bit pattern b; is given by 
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X, (bi) = E{X; [bi} 
M 
= = A7TbR + 1S clare 4 +2 YCEA*TH; 
k=] 


4int 


(19) 


From equation (18) and equation (BS) of Appendix B, the conditional variance ox (bi) 
of X; can be calculated as follows: 


of, (Bi) = S.CRA*T0% (D2. +08, o}frv/4x?(0? +82) 
k=l 


~Gashparsagp OP lea 
—e2™ (2nvsin 2nd, + 2x5, cos2n5, I (20) 
+4n?(85 -v’ J2nv -4x2(v? + 5) 

+e?" (2x5, sin2n5, - 2xvcos2n5,)}} 


where v=BT and 5, =@,T/2x . Here f is the laser linewidth. In summary, the 


Gaussian approximation allows the decision variable Y; in equation (2) to be considered 
as a Gaussian random variable with conditional mean Yi(bi) = Xj bi) in equation (19) 


and conditional variance OY, (bi) given by 


o¥, (bi) = ox, bi) + ow (bi)+ ow, + ON. (21) 


For a threshold a, the conditional bit error probability given bit patterns bY = (doo, 
=1 % 
bi1,-15°"» Dim,-1, bit,0,°-,bim,o) and bj = (O10, bi,-1.°--, Ding,-1 Bj1,0.°~im.o) is [Ref. 20] 


Ph (Dit,-1>---»Bing,-1>Bi2,0r---»bigg 0) = =o(bt) + =h(ti) (22a) 
where 
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where 


sae a-¥ji (b?) 
Po(6: | = Serf in, =) (22b) 
m1\_ 1 ¥;(bi)- a 
P,(Bi) = scr Fao, (22c) 
and erfc( - ) is again defined (as in equation (10d)) as 
erfo(a) = = [re ax. (224) 


The bit error probability P, is obtained by taking the expectation of P,(bj1_-1,--, bim.-1, 
bi1,0,°--.biM,0) With respect to the bit patterns (b;),.1,°*-, Dim,-15 Di1,0."":Dim.0)- Since there are 


22M such patterns, the bit error probability becomes 


1 
PB, = FEM Lieb ctenbanes .bi3.0 bao) b(Dit,-1>-» Ding ~1>Di1,09--> bim,o) (23) 


gooey 


where the summation is over all 22M patterns (bj;,-1,°--, bim,-1, bi1,0°°°,bim,o). The optimal 
threshold that minimizes the bit error probability is the value that satisfies equation (24). 
= (-0\? - 
1_-[a-¥i(82)} /203, (62) 


2 +1 ot gece Dang apo Diy Queer Dj -0 
(D1, -1 seo Dig —1>Di1,0----Dimo) oy, (7) 


(24) 


7 are l® (af /20%, (8) |=, 


Equation (24) is obtained by setting P,/da. to zero. 
In the case when all laser sources have the same mean wavelength but are 


uncorrelated, the above results apply by setting 5, =0 , k = 1, 2,---,M in equation (21). 
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Furthermore, when all laser sources are identical such that all sources have the same 


wavelength and phases noise process except for the random initial phases, then OX. (bi) 


in equation (20) reduces to 


M 
o% (bi) = 2 LCA Toh +b%o). (25) 


B. NUMERICAL ANALYSIS 

In this section, numerical results are presented for a syste with a bit rate of 500 
Mb/s. The responsitivity of the photodiode is taken to be 0.5 and the laser extinction 
ratio is 1/20. The dark current Iy, = 10 nA. The effective noise temperature of a low 
noise amplfier-integrate and dump-slicer receiver is 180 K. Furthermore, assuming a 
matched load of Ry = 50 Q, the post-detection thermal noise current spectral density is 
No = 2kTo/R, = 1022 A?/Hz [Ref. 15 and Ref. 20], where k = 1.38 x 10-23 J/K is the 
Boltzmann's constant. 

This section is broken down into the three cases for the probability of bit error plots 
and the power penalty plot summaries. The power penalty plots compare the case of two 
adjacent channels having synchronous transmission with the case of two adjacent 
channels having asynchronous transmissicn. The cases summarized in the three power 
penalty plots are for a bit error probability Pp of 10-15. The data collected is the result of 
the MATLAB model of the optical chip interconnect for asynchronous transmission 
(Appendix D) and synchronous transmission (Appendix C). The first case of laser source 
occurs when all channels operate with independent laser sources. Case two occurs when 
all laser sources have the same mean wavelength but have different noise processes, and 
case three occurs when all laser sources are identical with the exception of tie initial 


phases. 
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1. Casel 

Figures 20-21 show the bit: »robability P,, versus the received peak power 
A?/2 for various levels of crosstalk from two adjacent channels relative to that of a single 
channel operation (zero crosstalk) and with laser linewidth-bit rate ratio and frequency 
spacing-bit rate ratios (v,5,,52) taken to be (0.1,0.3,0.3) and (0.1,0.7,0.7), respectively. 
From Fig. 20, the bit error rate floor exists around P, = 10°!5 for—20 dB crosstalk 
irrespective of the received peak power. Crosstalk levels must be less than —26 dB for a 
power penalty of 1 dB or less for P, = 10°15, From Fig. 21, it is seen that by increasing 


the frequency the frequency spacing the crosstalk can be reduced to less than —23 dB for 
1 dB or less in power penalty at P, = 10°15. 
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Figure 20: P) vs. A2/2 as a function of crosstalk levels with (v,5,,52) = (0.1,0.3,0.3). 
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Figure 22: Normalized optimal threshold a vs. A2/2 as a function of crosstalk levels 
with (v,5,,52) = (0.1, 0.3, 0.3). 
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Figure 23: Normalized optimal threshold a vs. A2/2 as a function of crosstalk levels 
with (v,53 552) = (0.1, 0.7, 0.7). 


Figures 22-23 show the normalized optimum threshold versus the received 
peak power as function of crosstalk levels. The optimum threshold decreases with 
increasing crosstalk. Figure 24 shows the results for (v,8,,55) = (1,0.3,0.3). It is seen 
that a laser with a larger laser linewidth-bit rate ratio improves the performance. For 0.7 
dB or less in power penalty when P, = 10-15, the permitted crosstalk level is less than 
—23 dB instead of —26 dB as in Fig. 20. When the normalized linewidth v is increased 
to 5 as in Fig. 25, there is less than 0.7 dB power penalty for crosstalk levels less then 
-—17 dB. This happens because only a portion of the crosstalk energy falls within the 
detection bandwidth of 500 MHz. This result encourages the use of lasers with a large v 
as long as the waveguide bandwidth is larger than the signal spectrum. When the signal 
spectrum broadened by the laser phase noise approaches the waveguide bandwidth, loss 


in signal begins to occur and performance deteriorates rapidly. 
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Figure 25: P vs. A2/2 as a function of crosstalk levels with (v,5,,52) = (5,0.3,0.3). 
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2. Case II 
Figure 26 shows the performance of case two when all laser sources have the 
same mean wavelength but with different noise process for (v,5,,52) = (0.1, 0, 0). A bit 
error rate floor exists at P,, = 10-14 irrespective of the received peak power. The permitted 
crosstalk level is less than —26 dB for a power penalty of 1 dB or less at P, = 10°15. In 
general, the performance of case two is always worse than that of case one given the 


same and crosstalk level. It is obvious that the crosstalk effect is reduced by using laser 


sources with different wavelengths. 
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Figure 26: Pp vs. A2/2 as a function of crosstalk levels with (v,5,,52) = (0.1,0,0). 
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3. Case Ill 
Figure 27 shows the performance of case three when all lasers are identical 
with the exception of the initial phases. The performance is slightly better than that in 
Fig. 26. These results relate directly back to the derivation of equation (25) and are 


opposite to those in Case III for synchronous transmission. 
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Figure 27: P, vs. A2/2 as a function of crosstalk levels with all laser sources having 
identical phase noise processes except for the initial phases. 


4. Power Penalty Plot Summary 
Figure 28 shows the power penalty versus crosstalk level for two interfering 
channels with a normalized frequency spacing of 0.3 as a function of the normalized laser 
linewidth. The power penalty for the synchronous transmission where 5, = 0.3 in 


equation (1) is also plotted for comparison. In general, the power penalty is smaller for 
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asynchronous transmission when the linewidth is less than 3 for crosstalk levels up to 
—20 dB. For larger laser linewidths, the power penalty for synchronous transmission is 
slightly less. Figure 29 shows the power penalty versus normalized frequency spacing 
for two interfering channels with a normalized linewidth of 0.1 as a function of the 
crosstalk level. In general, the power penalty for asynchronous transmission is less than 
that of synchronous transmission. Figure 30 shows the power penalty for two interfering 
channels versus the normalized laser linewidth with a normalized frequency spacing of 
0.3, as a function of the crosstalk level. In general, the power penalty for asynchronous 
transmission for normalized laser linewidths less than 3 is less than that for synchronous 


transmission. 





Figure 28: Power Penalty vs. crosstalk level for a normalized frequency spacing of 
5, = 5, = 0.3 for asynchronous and synchronous transmission. 
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Figure 29: Power penalty vs. normalized frequency spacing (5, = 52) for a 
normalized linewidth of v = 0.1 for asynchronous and synchronous transmission. 


Power Penalty (dB) 








Figure 30: Power penalty vs. normalized linewidth for 5, = 5, = 0.3 
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C. DISCUSSION 

To conclude, this chapter presented a mathematical framework to analyze the 
performance of asynchronous optical chip interconnects in terms of the bit error 
probability versus the received peak power as a function of the crosstalk level, frequency 
spacing, and laser linewidth. Furthermore, the analysis can handle any number of 
adjacent channels. The conclusion drawn from this investigation is that adjacent channels 
must use laser sources of different wavelengths to reduce the effect of crosstalk. Laser 
sources with a large linewidth also help, as long as the waveguide bandwidth is much 
larger than the resulting signal spectrum. In fact, this is the only way to reduce the effect 
of a given crosstalk level when all laser sources are locked to a master source. The 
performance depends explicitly on the crosstalk level given a received peak power. 
Finally, as seen by Figs. 28, 29, and 30, asynchronous optical chip interconnects perform 


better than their synchronous counterparts. 

















VI. CONCLUSION 


This thesis derived a crosstalk model for an optical chip interconnect that determines 
the crosstalk-induced system penalty caused by the coupling between two or more 
adjacent channels in a single-mode waveguide array. MATLAB version 3.5k was useful 
in the development and the verification of the model. The resulting performance curves 
and power penalty plot summaries demonstrate the effectiveness of the model for both 
synchronous and asynchronous transmission. In general, the asynchronous case proved 
to be more sensitive than the synchronous case for crosstalk levels less than —20 dB and 
normalized laser linewidths less than 2. These results follow those in Ref. 9 where it was 
determined that the sensitivity of an asynchronous receiver is only about 0.5 dB less than 
its synchronous counterpart. 

Study is currently underway to determine the effect on this crosstalk model of 
shutting down one or more of the adjacent channels during transmission. The results 


from this study should further verify the effectiveness of this model. 
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APPENDIX A - DERIVATION OF THE CONDITIONAL 
MEAN SQUARE VALUE OF X, GIVEN b; FOR 
SYNCHRONOUS TRANSMISSION 


The conditional mean square value of X; given bi is important in the derivation of the 
conditional variance o%, (b i) of X; (equation (8)). The derivation of the conditional mean 


square begins by finding the expectation of X? (see equation (6)) given bi. 
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In the derivation of equation (A.1), all the initial phase differences, >, k = 1, 2,°--, 
M are uniform variables over (0, 27). 

Let y= Oo(t) — 0, (t) —[@9(t) — ©, (t)]. The laser phase noise 6,(t), k = 0, 1,---, M 
is characterized by a Wiener process [Ref. 17 - Ref. 18] such that dO, (t)/dt = 2my(t) 
where ji(t) is a zero mean white Gaussian process of PSD B/2x Hz where B is the laser 
linewidth (assumed to be the same for all laser sources). The variance of 0,(t) is 2Bt. 


The process y is a zero mean Gaussian process and can be expressed as 


y= 2nf [Holt )- H(t) By - 2nf [uo(t) — H(t) Bt 
= 2nJ [Ho(t) = Hy (ty)}aty . 


(A.2) 
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Therefore, the variance 0 is given by 
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Oy 
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= anf Po(, bod ty }dtydt, 


4nB fi (ita -lo)8(u)du 
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(A.3) 


Using the fact that y and —w are both Gaussian random variable with zero mean and 


variance 4nB|t — 1], the expectation value is [Ref. 13] 


E{ei¥} = E{e¥} = CWP? xg 2nBit tl 


(A.4) 


Substituting equation (A.4) into equation (A.1), the conditional mean variance is simplified 
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where v= BT and 6, = @,T/2n. 
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APPENDIX B - DERIVATION OF THE CONDITIONAL MEAN 
SQUARE VALUE OF X, GIVEN b; FOR ASYNCHRONOUS 
TRANSMISSION 


The conditional mean square value of X; given bj is important in the derivation of 
the conditional variance o (bi) of X; (equation (20)). The random variable X; is 


derived from equations (14) and (15b). By taking the expectation of X? given bi and t 
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The conditional variance of X; given by bj and 7 is given by 


o%, (5.2) 
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It should be noted that the derivation of equation (B.4) depends on the fact that all initial 
phase differences , with k = 1, 2,...,M are uniform random variables over (0,27). 

Let y = Oo(t)— ©, (t)—[09(t)— 0, (t)]. The laser phase noise @,(t) with k = 0, 1,--, 
M is characterized by a Wiener process [Ref. 17 - Ref. 18] such that d@,(t)/dt = 2x1(t) 
where 1(t) is a zero mean white Gaussian process of PSD B/2x Hz where f is the laser 
linewidth (assumed to be the same for all laser sources). The variance of 0,(t) is 2nBt. 


The process y is a zero mean Gaussian process and can be expressed as 
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y= 2nf [Ho(u)- H(t) At - 2nf [uo(t) — H(t, Bh 
= 2x} [Ho(t:)- Hx (th)]dt, . @.5) 


Therefore, the variance 0%, is given by 
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Using the fact that y and —y are both Gaussian random variables with zero mean and 


variance 4nxfjt—~ 1], the expected value is [Ref. 13] 


Substituting equation (B.7) into (B.4), the conditional variance obtained is 
2 b; p= =a A‘ M C2 2 b2 Tt, rt —2nfi1—7 dtd 
0x, (bi,t) = > Lerbiobi.f, "Joe cosa, (t~ t)dtdt 
k=1 
M T (T 
+2 Chbiobicol, [een cos, (t — swt} 
k=1 
A‘ | M . 
Al Scheid. (ty —lul)e?*""! coso,udu (B.8) 
k=l 


M = 
+). Cibjobi..o hee ((T- 2, |-lu))e 7" cos, vs} 
k=1 


Finally, the conditional variance of X; given bi is obtained by averaging Of, (bi,7) 
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over T giving the end result as 


- M 
oi,(B) = Detatren(eha+bhe Nao ail 
1 
64x° (v? + 53) 
x(8x7v5, [2n8, —e°* (2nvsin 2nb, +25, cos2n5,)| 
+4n?(83 - v?) (B.9) 
anv - 4nx?(v? + 8) +e" 


x(2n5, sin2x5, — 2nvcos275, )))} 


where v= BT and 5, = @,T/2n. 
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APPENDIX C - MATLAB MODEL OF OPTICAL CHIP 
INTERCONNECT FOR SYNCHRONOUS TRANSMISSION 


This model is for the case of 2 adjacent channels with (v,5,,52) = (0.1,0.3,0.3) for 
synchronous transmission. The MATLAB version used was version 3.5k. The program 
itself is a translation of equations (1) through (13) into MATLAB code. 


% Synchronous Case: No=1E-22 

% Case 1.3: nu=.1 and delta function subscript k (delk)=.3 
% th2ch13x.m 

% 19 July 93 

% Known constants 


R=.5; 

T=1/(500E6); 

NO=1E-22; 

q=1.6E-19; 

Idk=10E-9; % dark current 


% Coupling constants for Ck*2 
ck1=0; 
ck2=3.16E-4; %-35 dB 


ck5=2.52E-3; %-26 


ck7=1 995E-2; %-17 dB 

ck8=3.981E-2; %-14 dB 

ck9=7.943E-2; %-11 dB 

ck10=.1585; %- 8 dB 

ck=[0 3.16E-4 6.32E-4 1.26E-3 2.52E-3 5.012E-3 .01 1.995E-2 3.981E-2 7.943E-2 
1585); 


% Pattern bit values 

%b1=[.2182 .9759 .9759 .2182]; %for b11,b01 for patterns A,B,C,D 
Gb2=[.9759 .2182 .9759 .2182}; %for b12,b02 for patterns A,B,C,D 
b1A=.2182; 

b1B=.9759; 

b1C=.9759; 

biD=.2182; 

b2A=.9759; 

b2B=.2182; 

b2C=.9759; 

b2D=.2182; 

b3A=b1/ *2+b2A42; 

b3B=b1b: :2+b2B“2; 

b3C=b1C42+b2C42; 

b3D=b1D42+b2D“2; 
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b00=.2182; % approximate zero 
b10=.9759; % approximate 1 


% ce acim k constant for the variance of X 

nu= 

delk=.3; 

ka=(2*pi*nu)/(4*pi*2*(nuA2+delk’2)); 

kb=1/(16*pi44*(nu42+delk42)2); 

aT Har re add bik iron etal tect A alone ¢ rar ht 
kd=(2*pi*delk*exp(-2*pi*nu))*(2*pi* (delk*cos(2*pi*delk)+nu*sin(2*pi*delk))); 

ke=4*piA2*(nuA2-delkA2); 

k=ka-kb*(kc+kd+ke); 

n=]; 

j=l, 

m=]; 


% -60dB<=Peak Power<=-44dB 
% Peak power is A%2/2 = Ps 
for ss=-60:.2:-44; 
ps(m)=10(ss/10); 
m=m-+1; 
end 


% The goal here is to solve the equation for the optimal threshold 
% vs peak power 


alpha=zeros(length(ps),length(ck)); 
PBE=zeros(length(ps),length(ck)); 





for j=1:length(ck); % coupling values loop 


n=]; 
for n=1:length(ps); % peak ppower values 
Ps=ps(n); 

% Constants 


qcnst=q/(Ps*T*R); 


vnoise=NO/(R“2*Ps*2*T); 
vdark=(q*Idk)/(R“2*Ps“*2*T); 


% Determine the sigma for each of the 0 bit patterns 
sOA=sqrt(4*ck(j)*b00*2*k*b3A+qcnst*(b00“2+ck(j)*b3A)+vnoise+vdark); 
sOB=sqrt(4*ck(j)*b00*2*k*b3B+qcnst*(b00*2+ck(j)*b3B)+vnoise+vdark); 
sO0C=sqrt(4*ck(j)*b0042*k*b3C+qenst*(b0042+ck(j)*b3C)+vnoise+vdark); 
sOD=sqrt(4*ck(j)*b00*2*k*b3D+qenst*(b00“2+ck(j)*b3D)+vnoise+vdark); 
% Determine the variance for each of the 0 bit patterns 


s0A2=s0A‘2; 
s0B2=s0B‘2; 
s0C2=s0C‘2; 
s0D2=s0D“2; 


% Determine the sigma for each of the 1 bit patterns 
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SS ee 


s1A=sqrt(4*ck(j)*b10*2*k*b3A +qcnst*(b102+ck(j)*b3A + vnoise+vdark); 
$1B=sqrt(4*ck(j)*b10°2*k*b3B+qcnst*(b10°2+ck(j)*b3B)+vnoise+vdark); 
$1C=sqri(4*ck(j)*b10°2*k*b3C-+qcnst*(b1042+ck(j)*b3C)+vnoise+vdark); 
$1D=sqrt(4*ck(j)*b10°2*k*b3D+qenst*(b10°2+ck(j)*b3D)+vnoise+vdark); 





% Determine the variance for each of the 1 bit patterns 


$1A2=s1A‘2; 
s1B2=s1 B42; 
s1C2=s1C*2; 
s1D2=s1D“2; 


% Determine the mean for each of the 0 bit patterns 


m0A=b0042+ck(j)*b3A; 
m0B=b002+ck(j)*b3B; 
m0C=b00A2+ck(j)*b3C; 
m0D=b002+ck(j)*b3D; 


% Determine the mean for each of the 1 bit patterns 


m1A=b1042+ck(j)*b3A; 
m1B=b1042+ck(j)*b3B; 
m1C=b10°2+ck(j)*b3C; 
m1D=b1042+ck(j)*b3D; 


%Solve the equation F(a)=0 to optimize a, the threshold. F(a)=0 comes 
% from (dP/da)=0 


a=0:.001:1; 
FA=(1/sOA)*exp(-((a-m0A).42)./(2*sOA2))-(1/s1A)*exp(-((m1A-a).42)/(2*s1A2)); 
FB=(1/sOB)*exp(-((a-m0B).42)./(2*s0B2))-(1/s1B)*exp(-((m1B-a).42)/(2*s1B2)); 
FC=(1/sO0C)*exp(-((a-m0C).42)./(2*s0C2))-(1/s1C)*exp(-((m1C-a).42)/(2*s1C2)); 
FD=(1/sOD)*exp(-((a-m0D).42)./(2*sOD2))-(1/s1D)*exp(-((m1D-a).2)/(2*s1D2)); 
F=FA+FB+FC+FD; 


[Y,=min(abs(F)); 


alpha(n,j)=a(1); 

aa=alpha(n,j); 

% Determine the probability of bit error for each optimal threshold 
% over all bit patterns 4,R C.D 
pOA(n)=erf((alpha(n,j)-m0A)/sOA nf); 
pOB(n)=erf((alpha(n,})-m0B)/sOB sinf); 
pOC(n)=erf((alpha(n,j)-m0C)/sOC, inf); 
pOD(n)=erf((alpha(n,j)-m0D)/sOD, inf); 


p1A(n)=erf((m1A-alpha(n,j))/s1A, inf); 
p1B(n)=erf((m1B-alpha(n,}))/s1B,inf); 


p1C(n)=erf((m1C-alpha(n,j))/s1C,inf); 
p1D(n)=erf((m1D-alpha(n,j))/s1D, inf); 


PBE(n,j)=(p0A(n)+p1 A (n)+p0B(n)+p1B(n)+p0C(n)+p1C(n)+p0D(n)+p1D(n))/16; 
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% n=n+1; 


end % Ps loop 
% j=j+1; 
end %Ck loop 


al=alpha(:,1); 
a2=alpha(:,2); 
a3=alpha(:,3); 
ad=alpha(:,4); 
a5=alpha(:,5); 
a6=alpha(:,6); 
a7=alpha(:,7); 
a8=alpha(:,8); 
a9=alpha(:,9); 
a10=alpha( 10); 
al l=alpha(:,11); 


p1=(PBE(:,1)); 
p2=(PBE(:,2)); 
E(:,3)); 


5=(PBE(: 5); 
p6=(PBEC 6)); 
p7=(PBE(:,7)); 
p8=(PBEC 8)); 
po=(PBE(.. 9)); 
p10=PBE(:,10); 
pll=PBE(:,11); 


PB=10*log10(ps); 
plot(PB,al,’-',PB,a2,’-',PB,a3,'-',PB,a4,'-',PB,a5,'-',PB,a6,'-',PB,a7,'-'),.. 


grid... 
%title(‘case 1.3 (nu=.1,delk=.3)’),.. 
xlabel(‘Peak Power (dBW)'),ylabel(Optimal Threshold’) 


s=['Crosstalk=']; 
text(.15,.6,s,'sc') 


c1=0, 

s=(' ‘num2str(c1)]; 
text(.2,.55,s,'sc' 
[xs,ys}=dc2sc(PB(48),a1(48)); 
polyline([.3,xs},[.56,ys],'-r','sc’) 


c2=-35; 

s=(num2str(c2),' dB']; 
text(.2,.52,s,'sc’) 
[xs,ys]=dc2sc(PB(48),a2(48)), 
polyline([.3,xs]},[.53,ys],'-1','sc’) 
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c3=-32; 

s=[num2str(c3),' dB’); 
text(.2,.49,s,'sc’ 
(xs,ys]=dc2sc(PB(48),a3(48)); 
polyline([.3,xs],[.5,ys],-r’,'sc’) 


c4z=-29; 

s=(num2str(c4),’ dB’; 
text(.2,.46,s,'sc’) 
{xs,ys}=dc2sc(PB(48),a4(48)); 
polyline([.3,xs},[.47,ys],'-r’,'sc’) 


c5=-26; 

s=[num2str(c5),' dB’); 
text(.2,.43,s,'sc’ 
[xs,ys}=dc2sc(PB(48),a5(48)); 
polyline((.3,xs),[.44,ys],'-r’,'sc’) 


c6=-23; 

s=(num2str(c6),' dB’); 
text(.2,.4,s,'sc’) 
[xs,ys]=dc2sc(PB(48),a6(48)); 
polyline({.3,xs),[.41,ys],'-r’,'sc’) 


c7=-20; 
erinnzenc). dB’); 
text(.2,.37,5, 
Ixs,ys}ede2se(PB(48),a7(48)); 
polyline([.3,xs],[.38,ys),'-1’,'sc’) 


meta t2ch13xa 


pause 


axis([-60 -44 -21 -5)) 

plot(PB,p1,’ -',PB,p2,'-',PB,p3,'-’,PB,p4,-',PB,p5,'-',PB,p6,'-',PB,p7,-’),.. 
semilogy,grid, 

Qtitle(Performance Curves: case 1.3 (nu=.1,delk=.3)’),.. 

xlabel(‘Peak Power (dBW)'),ylabel(Probability of Bit Error) 


s=['Crosstalk=]; 
text(.15,.6,s,'sc’) 


cl=-20; 

s={[num2str(c1),' dB’); 
text(.2,.55,s,'sc’) 
[xs,ys]=dc2sc(PB(38),p7(38)); 
polyline([.3,xs],[.56,ys],'-1’,'sc’) 


c2=-23; 
s=[num2str(c2), dB’); 
text(.2,.52,,'sc’) 





[xs,ys}=dc2sc(PB(38),p6(38)); 
polyline([.3,xs),{.53,ys),'-r’,'sc’) 


c3=-26; 

s=(num2str(c3),' dB’); 
text(.2,.49,s,'sc’) 
(xs,ys]=de2sc(PB(38),p5(38)); 
polyline({.3,xs),{.5,ys],'-1,'sc’) 


cA=-29; 

s=(num2str(c4),’ dB’); 
text(.2,.46,s,'sc') 
(xs,ys]}=dc2sc(PB(38),p4(38)), 
polyline({.3,xs),[.47,ys],'-1’,'sc’) 


c5=-32; 

s=[num2str(c5),' dB’); 
text(.2,.43,s,'sc’) 
(xs,ys]=dce2sc(PB(38),p3(38)); 
polyline([.3,xs),[.44,ys],'-1','sc’) 


c6=-35; 

s=[num2str(c6),' dB’); 
text(.2,.4,s,'sc’) 
[xs,ys]=dce2sc(PB(38),p2(38)); 
polyline([.3,xs],[.41,ys],'-1’,'sc’) 


c7=0; 

s=[' 'num2st(c7)]; 
text(.2,.37,s,'sc’) 
(xs,ys]=dc2sc(PB(38),p1(38)); 
polyline({.3,xs],[.38,ys],-1','sc’) 
axis; 


meta t2ch13xb 
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APPENDIX D - MATLAB MODEL OF AN OPTICAL CHIP 
INTERCONNECT FOR ASYNCHRONOUS TRANSMISSION 


This model is for the case of 2 adjacent channels with (v,5,,5,) = (0.1,0.3,0.3) for 
asynchronous transmission. The MATLAB version used was version 3.5k. The program 
itself is a translation of equations (14) through (25) into MATLAB code. 


% Asynchronous Case: No=1E-22 

% k=2 implying two channels 

% Case 1.3: nu=.1 and delta function subscript k (delk)=.3 
% th2as13x.m 

% 20 July 93 


% Known constants 


R=.5; 
T=1/(500E6); 
NO=1E-22; 
q=1.6E-19; 
Idk=10E-9; 


% Coupling constants for CkA2 
ck1=0; 

ck2=3.16E-4; %-35 dB 

ck3=6.32E-4; %-32 dB 

ck4=1.26E-3; %-29 dB 

ck5=2.52E-3; %-26 dB 

ck5a=5.011E-3; %-23 dB 

ck6=.01; %-20dB 

ck7=1.995E-2; %-17 dB 

ck8=3.981E-2; %-14 dB 

ck9=7.943E-2; %-11 dB 

ck10=.1585; %-8 dB 

sort ea 6.32E-4 1.26E-3 2.52E-3 5.011E-3 .01 1.995E-2 3.981E-2 7.943E-2 
el : 


% Pattern bit values 


%b1=[.2182 .9759 .9759 .2182 .2182 .2182 .9759 .9759 .2182 .2182 .2182 .9759 
9759 9759 .2182 .9759); 

%b2=[.9759 .2182 .9759 .2182 .2182 .9759 .2182 .9759 .2182 .2182 .9759 .2182 
9759 2182 .9759 .9759]); 


%for b32,b22,b12,b02 for patterns A,B,C,D,E,F,G,H,I,J,K,L,M,N,O,P 


%b3=[.9759 .9759 .9759 .9759 .2182 .2182 .2182 .2182 .2182 .9759 .2182 .2182 
2182 .9759 .9759 .9759]; 

%for b33,b23,b13,b03 

Fb4=[.9759 9759 .9759 .9759 .9759 .9759 .9759 .9759 .2182 .2182 .2182 .2182. 
2182 .2182 .2182 .2182); 
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%for b34,b24,b14,b04 


blA=.2182; 
b1B=.9759; 
b1C=.9759; 
b1D=.2182; 
blE=.2182; 
bl F=.2182; 

b1G=z.9759, 
blH=.9759; 
blI=.2182; 

blJ=.2182; 

b1K=.2182; 
b1L=.9759; 
b1M=.9759; 
b1N=.9759; 
b10=.2182; 
b1P=.9759; 

b2A=.9759; 
b2B=.2182; 
b2C=.9759; 
b2D=.2182; 
b2E=.2182; 
b2F=.9759; 

b2G=.2182; 
b2H=.9759; 
b2I=.2182; 

b2J=.2182; 

b2K=.9759; 
b2M=.9759; 
b2N=.2182; 
b20=.9759; 
b2P=.9759; 

b3A=.9759; 
b3B=.9759; 
b3C=.9759; 
b3D=.9759; 
b3E=.2182; 

b3F=.2182; 

b3G=.2182; 
b3H=.2182; 
b3I=.2182; 

b3J=.9759, 

b3K=.2182; 
b3L=.2182; 
b3M=.2182; 
b3N=.9759; 
b30=.9759; 
b3P=.9759; 

b4A=.9759; 
b4B=.9759; 
b4C=.9759; 
b4D=.9759; 
b4E=.9759; 
b4F=.9759; 
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b4G=.9759; 
b4H=.9759; 
b4I=.2182; 

b4J=.2182; 

b4K=.2182; 
c4+L=.2182; 
o4M=.2182; 
h4N=.2182; 
b40=.2182; 
b4P=.2182; 


bSA=b1A42+b2A42+b3A42+b4A‘2; 
b5B=b1B42+b2B42+b3B“2+b4B2; 
b5C=b1CA2+b2C42+b3C42+b4C%2; 
b5D=b1D42+b2D‘2+b3D42+b4D 2; 
bSE=b1E42+b2E42+b3E42+b4E“2; 
b5F=b1F42+b2F42+b3F42+b4F*2; 
b5G=b1G42+b2G42+b3G%2+b4G“2; 
bSH=b1H‘2+b2H‘42+b3H‘42+b4H“2; 
b5Il=b1H42+b2H‘2+b3H“2+b4142, 
b5J=b1J42+b2JA2+b3)42+b4542; 
bSK=b1K42+b2K42+b3K‘42+b4K‘2; 
b5M=b1M42+b2M42+b3M42+b4M“2, 
b5SN=b1N42+b2N42+b3N‘42+b4N“2; 
b5O=b1042+b2042+b3042+b40%2; 
b5P=b1P42+b2P42+b3P42+b4P42; 
b0O=.. . 52; % approximate zero 
b10=.9759; % approximate 1 


nu=.1; 
aoe 3; 
a=(pi*nu)/(4*pit2*(nuA2-+delk’2)); 
Kees" 146*(nuA2+delk42)43); 
ke=2*pitnutsin(2*pitdelk)+3¥pitdelk*cox(2ep*delk 
d=(8*pi42*nu*delk)*(2*pi*delk-(exp(-2*pi*nu))*kc); 
Kent" biedelky*sin(*yi*aelk). (2*pi*nu)*cos(2*pi*delk); 
kf=4*pi42*(delk*2-nu’2)*(2*pi*nu-4*pi*2*(nu42+delk‘2)+exp(-2*pi*nu)*ke); 
k=ka-kb*(kd+kf); 
n=1; 
j=l, 
m=1; 


% -60dB<=Peak Power<=-44dB 
% Peak power is A42/2 = Ps 
for ss=-60:.5:-44; 
ps(m)=10(ss/10); 
m=m+1; 
end 


% The goal here is to solve the equation for the optimal threshold 
% vs peak power 


alpha=zeros(length(ps),length(ck)); 








PBE=zeros(length(ps),length(ck)); 





for Senge: % coupling values loop 


n=1, 

for n=1:length(ps); % peak power values 
Ps=ps(n); 

% Constants 


qcnst=q/(Ps*T*R); 
vnoise=NO/(R4A2*Ps42*T); 
vdark=(q*Idk)/(R42*Ps42*T); 


% Determine the sigma for each of the 0 bit patterns 


sOA=sqrt(4*ck(j)*b0042*k*b5A +qcnst*(b0042+ck(j)/2*b5A)+vnoise+vdark); 
sOB=sqrt(4*ck(j)*b00“2*k*b5B+qcnst*(b00“2+ck(j)/2*b5B)+vnoise+vdark); 
sOC=sqrt(4*ck(j)*b00A2*k*b5C+qcnst*(b00*?+ck(j)/2*b5C)+vnoise+vdark); 
sOD=sqrt(4*ck(j)*b00*2*k*b5D+qcn st*(b0042+ck(j)/2*b5D)+vnoise+vdark); 
sOE=sqrt(4*ck(j)*b0042*k*bS5SE+qcnst*(b00‘2+ck(j)/2*b5E)+vnoise+vdark); 
sOF=sqrt(4*ck(j)*b00“2*k*b5F+qcnst*(b00‘2+ck(j)/2*bS5F)+vnoise+vdark); 
sOG=sqrt(4*ck(j)*b00*2*k*b5G+qcnst*(b00“2+ck(j)/2*b5G)+vnoise+vdark); 
sOH=sqrt(4*ck(j)*b0042*k*b5H+qcnst*(b002+ck(j)/2*bSH)+vnoise+vdark); 
sOI=sqrt(4*ck(j)*b00*2*k*bS5I+qcnst*(b0042+ck(j)/2*b51)+vnoise+vdark), 
sOJ=sqrt(4*ck(j)*b0042*k*b5J+qcnst*(b00“2+ck(j)/2*b5J)+vnoise+vdark); 
sOK=sqrt(4*ck(j)*b00“2*k*b5K+qcnst*(b00‘2+ck(j)/2*b5K)+vnoise+vdark); 
sOL=sqrt(4*ck(j)*b00°2*k*b5L+qcenst*(b002+ck(j)/2*b5L)+vnoise+vdark); 
sOM=sqrt(4*ck(j)*b00“2*k*b5M+qcnst*(b00%2+ck(j)/2*b5M)+vnoise+vdark); 
sON=sqrt(4*ck(j)*b00*2*k*b5N+qcnst*(b0042+ck(j)/2*b5N)+vnoise+vdark); 
s0O=sqrt(4*ck(j)*b00°2*k*b5O+qcenst*(b00%2+ck(j)/2*b50)+vnoise+vdark); 
sOP=sqrt(4*ck(j)*b002*k*b5P+qcenst*(b00*2+ck(j)/2*b5P)+vnoise+vdark); 


% Determine the variance for each of the 0 bit patterns 


s0OA2=s0A‘2; 
s0OB2=s0B“2; 
s0C2=s0C‘2; 
s0OD2=s0D‘2; 
sOE2=s0E“2; 
sOF2=sOF‘2; 
s0G2=s0G‘2; 
s0OH2=s0H“2; 
s012=s0I42; 
sOJ2=s0J42; 
sOK2=s0OK‘2; 
sOL2=s0L‘2; 
sOM2=s0M“2; 
sON2=sON‘“2; 
s0O2=s00“2; 
sOP2=sOP“2; 


% Determine the sigma for each of the 1 bit patterns 


s1A=sqrt(4*ck(j)*b10*2*k*b5A+qenst*(b1042+ck(j)/2*b5A)+vnoise+vdark); 
s1B=sqrt(4*ck(j)*b1042*k*b5B+qcnst*(b1042+ck(j)/2*b5B)+vnoise+vdark); 
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$1C=sqrt(4*ck(j)*b10°2*k*b5SC+qcnst*(b1042+ck(j)/2*b5C)+vnoise+vdark); 
$1D=sqrt(4*ck(j)*b10°2*k*bS5D+qcenst*(b102+ck(j/2*b5D)+vnoise+vdark); 
slE=sqrt(4*ck(j)*b10°2*k*bSE+qcnst*(b1042+ck(j)/2*bSE)+vnoise+vdark), 
$1F=sqrt(4*ck(j)*b10A2*k*bSF+qcnst*(b102+ck(})/2*b5F)+vnoise+vdark); 
s1Gzsqrt(4*ck(j)*b10°2*k*b5G-+qcnst*(b102+ck(j)/2*b5G)+vnoise+vdark); 
s1H=sqrt(4*ck(j)*b10°2*k*b5H+qcnst*(b10%2+ck(j)/2*b5H)+vnoise+vdark); 
sll=sqrt(4*ck(j)*b1042*k*bSI+qcnst*(b1042+ck(j)/2*b51)+vnoise+vdark); 

sl J=sqrt(4*ck(j)*b10°2*k*b5J+qcenst*(b1042+ck(j)/2*b5J)+vnoise+vdark); 
sl Kesqrt(4*ck(j)*b10*2*k*bSK+qcnst*(b10°2+ck(j)/2*b5K)+vnoise+vdark); 
slL=sqrt(4*ck(j)*b10*2*k*b5L+qcnst*(b10A2+ck(j/2*b5L)+vnoise+vdark); 
s1Me=sqrt(4*ck(j)*b1042*k*b5M+qcnst*(b10°2+ck(j)/2*b5M)+vnoise+vdark); 
s1N=sqrt(4*ck(j)*b10°2*k*b5N+qcnst*(b10%2+ck(j/2*b5N)+vnoise+vdark); 
s10=sqrt(4*ck(j)*b10°2*k*b5O+qcnst*(b102+ck(j)/2*b5O)+vnoise+vdark); 
$1P=sqrt(4*ck(j)*b10°2*k*b5P+qcnst*(b102+ck(j)/2*b5P)+vnoise+vdark); 


% Determine the variance for each of the 1 bit patterns 


s1A2=s1A42; 
$1B2=s1B42; 
s1C2=s1C*2; 
s1D2=s1D“2; 
slE2=slE‘2; 
s1F2=s1F42; 
$1G2=s1G42; 
s1H2=s1H‘2; 
s1I2=s1142; 
s1J2=s1J42; 
s1K2=s1K42; 
s1L2=s1L42; 
s1M2=s1M‘42; 
sIN2=s1N42; 
s102=s10/2; 
s1P2=s1P2; 


% Determine the mean for each of the 0 bit patterns 


m0A=b002+ck(j)/2*b5A; 
m0B=b0042+ck(j)/2*bSB; 
m0C=b002+ck(j)/2*b5C; 
m0D=b002+ck(j)/2*bS5D; 
m0E=b002+ck(j)/2*bSE; 
m0F=b0042+ck(j)/2*b5F; 
m0G=b0042+ck(j)/2*b5G; 
m0H=b00“2+ck(j)/2*b5H; 
m0I=b00“2+ck(j)/2*b5I; 
m0J=b00°2+ck(j)/2*b5J; 
m0K=b0042+ck(j)/2*b5SK; 
mOL=b00“2+ck(j)/2*b5L; 
m0OM=b0042+ck(j)/2*b5M; 
mON=b00%2+ck(j)/2*bSN; 
m0O0=b00“2+ck(j)/2*b50; 
m0OP=b0042+ck(j)/2*bSP; 
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% Determine the mean for each of the 1 bit patterns 


m1A=b10/2+ck(j)/2*b5A; 
m1B=b1042+ck(j)/2*b5B; 
m1C=b1042+ck(j)/2*b5C; 
m1D=b1042+ck(j)/2*bSD; 
mlE=b10A2+ck(j)/2*bS5E; 
m1F=xb102+ck(j)/2*b5F; 
m1G=b1042+ck(j)/2*b5G; 
m1H=b1042+ck(j)/2*b5H; 
m1]=b1042+ck(j)/2*b5I; 
m1J=b1042+ck(j)/2*bS5J; 
m1K=b10/2+ck(j)/2*b5K; 
m1L=b1042+ck(j)/2*b5L; 
m1M=b1042+ck(j)/2*b5M; 
m1N=b1042+ck(j)/2*b5N; 
m10=b1042+ck(j)/2*b50O; 
m1P=b1042+ck(j)/2*b5P; 


%Solve the equation F(a)=0 to optimize a, the threshold. F(a)=0 comes 
% from (dP/da)=0 


a=0:.0009:1; 
FA=(1/sOA)*exp(-((a-m0A).‘2)./(2*sOA2))-(1/s1A)*exp(-((m1 Aa).42)/(2*s1A2)); 
FB=(1/sOB)*exp(-((a-m0B).42)./(2*s0B2))-(1/s1B)*exp(-((m1B-a).42)/(2*s1B2)); 
FC=(1/sOC)*exp(-((a-m0C).42)./(2*s0C2))-(1/s1C)*exp(-((m1C-a).42)/(2*s1C2)); 
=(1/sOD)*exp(-((a-m0D).42)./(2*s0D2))-(1/s1D)*exp(-((m1D-a).42)/(2*s1D2)); 
=(1/sOE)*exp(-((a-m0E).42)./(2*s0E2))-(1/s1E)*exp(-((m1E-a).42)/(2*s1E2)); 
FF=(1/sOF)*exp(-((a-m0F).‘2)./(2*sOF2))-(1/s1F)*exp(-((m1F-a).42)/(2*s1F2)); 
FG=(1/s0G)*exp(-((a-m0G).42)./(2*s0G2))-(1/s1G)*exp(-((m1G-a).42)/(2*s1G2)); 
FH=(1/sOH)*exp(-((a-m0H).42)./(2*sOH2))-(1/s1H)*exp(-((m1H-a).42)/(2*s1H2)); 
=(1/s0I)*exp(-((a-m0I).42)./(2*s012))-(1/s1I)*exp(-((m11-a).42)/(2*s112)); 
FJ=(1/s0J)*exp(-((a-m0J).2)./(2*s0J2))-(1/s1J)*exp(-((m1J-a).42)/(2*s1J2)); 
FK=(1/sOK)*exp(-((a-m0K).42)./(2*sOK2))-(1/s1K)*exp(-((m1K-a).42)/(2*s1K2)); 
FL=(1/sOL)*exp(-((a-m0L).42)./(2*s0L2))-(1/s1L)*exp(-((m1L-a).42)/(2*s1L2)); 
=(1/sOM)*exp(-((a-m0M).42)./(2*sOM2))-(1/s1M)*exp(-((m1M- 
a).42)/(2*s1M2)); 
FN=(1/sON)*exp(-((a-m0ON).42)./(2*sON2))-(1/s1N)*exp(-((m1N-a).42)/(2*s1N2)); 
FO=(1/s0O)*exp(-((a-m00O).42)./(2*s002))-(1/s10)*exp(-((m1O-a).42)/(2*s102)); 
FP=(1/sOP)*exp(-((a-m0P).*2)./(2*sOP2))-(1/s1P)*exp(-((m1P-a).42)/(2*s1P2)); 
F=FA+FB+FC+FD+FE+FF+FG+FH+FI+FJ+FK+FL+FM+FN+FO+FP; 


[Y,T}=min(abs(F)); 


alpha(n,j)=a(1); 
aa=alpha(n,j); 
% Determine the probability of bit error for each optimal threshold 
% over all bit patterns A,B,C,D 
pOA(n)=erf((alpha(n,j)-m0A)/sOA, inf); 
pOB(n)=erf((alpha(n,j)-m0B)/sOB, inf); 
p0C(n)=erf((alpha(n,j)-m0C)/sOC, inf); 
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pOD(n)=erf((alpha(nj)-m0D)/sOD, inf); 
pOE(n)=erf((alpha(n,j)-m0E)/sOE. inf); 
pOF(n)=erf((alpha(n,})-mOF)/sOF, inf); 
p0G(n)=erf((alpha(n,j)-m0G)/s0G, inf); 
pOH(n)=erf((alpha(n,j)-m0H)/sOH, inf); 
pOl(n)=erf((alpha(n,j)-m0I)/sOL inf), 
irr nb ns wrest ard cet 
pOK(n)=erf((alpha(n,j)-m0K)/sOK, inf); 
pOL.(n)=erf((alpha(nJ)-m0L)/sOL, inf); 
pOM(n)=erf((alpha(n,j)-m0M)/sOM, inf); 
pON(n)=erf((alpha(n,j)-mON)/sON, inf); 
p0O(n)=erf((alpha(n,j)-m0O)/s0O, inf); 
pOP(n)=erf((alpha(n,j)-m0P)/sOP, inf); 


p1A(n)=erf((m1A-alpha(n,j))/s1A,inf); 
p1B(n)=erf((m1B-alpha(n,}))/s1B, inf); 
p1C(n)=erf((m1C-alpha(n,))/s1C, inf); 
p1D(n)=erf((m1D-alpha(n,j))/s1D, inf); 
p1E(n)=erf((m1E-alpha(n,j))/s1E,inf); 
p1F(n)=erf((m1F-alpha(n))/s1F, inf); 
p1G(n)=erf((m1G-alpha(n,j))/s1G, inf); 
p1H(n)=erf((m1H-alpha(n,j))/s1H, inf); 
p1I(n)=erf((m11-alpha(n,j))/s1Linf); 
p1J(n)=erf((m1J-alpha(n,j))/s1J,inf); 
p1K(n)=erf((m1K-alpha(n,j))/s1K,inf); 
p1L(n)=erf((m1L-alpha(n,}))/s1L,inf); 
p1M(n)=erf((m1M-alpha(n,j))/s1M,inf); 
p1N(n)=erf((m1N-alpha(n,j))/s1N, inf); 
p10(n)=erf((m1O-alpha(n,j))/s10, inf); 
p1P(n)=erf((m1P-alpha(n,j))/s1P, inf); 


pA(n)=p0A(n)+p1A(n); 
pB(n)=p0B(n)+p1B(n); 
pC(n)=p0C(n)+p1C(n); 
pD(n)=pOD(n)+p1D(n); 
pE(n)=p0E(n)+p1E(n); 
pF(n)=p0F(n)+p1F(n); 
pG(n)=p0G(n)+p1G(n); 
pH(n)=p0H(n)+p1H(n); 
pl(n)=pOl(n)+p11(n); 
pJ(n)=p0J(n)+p1J(n), 
pK(n)=p0K(n)+p1K(n); 
pL{n)=p0L(n)+p1L(n); 
pM(n)=p0M(n)+p1M(n); 
pN(n)=pON(n)+p1N(n); 
pO(n)=p0O(n)+p10(n); 
pP(n)=p0P(n)+p1P(n), 


PBE(n,j)=(pA(n)+pB(n)+pC(n)+pD(n)+pE(n)+pF(n)+pG(n)+pH(n)+pI(n)+pJ(n)+pK 
(n)+pL(n)+pM(n)+pN(n)+pO(n)+pP(n))/64; 


end % Ps loop 








end %Ck loop 
al=alpha(:,1); 


2,6 
a7=alpha(:,7); 
a8=alpha(:,8); 
a9=alpha(:,9); 
alQ=alpha(:,10); 
al l=alpha(:,11); 


p2=(PBE(:,2)); 
p3=(PBE(:,3)); 


p7=(PBE(:,7)); 
p8=(PBE(:,8)); 
p9=(PBE(:,9)); 
p10=PBE(:,10); 
pll=PBE(:,11); 


PB=10*log10(ps); 
ae ae -',PB,a2,'-',P B,a3 Ae B a4,'-",P B,a5,'-',P. B ,26,-',P. B,a7,'~'),.. 


grid... 
%title(‘case 1.3,nu=.1,delk=.3: Four channels-Asynchronous’),.. 
xlabel(‘Peak Power(dBW)’),ylabel(‘Optimal Threshold’) 


s=['Crosstalk="']; 
text(.15,.6,s,'sc’) 


cl=0; 

s=[' ‘num2str(cl)); 
text(.2,.55,S,'sc’) 
[xs,ys}=dc2sc(PB(19),a1(19)); 
polyline(([.3,xs],[.56,ys],'-1','sc’) 


c2=-35; 

s=[num2str(c2),' dB’); 
text(.2,.52,s,'SC') 
[xs,ys]=dc2sc(PB(19),a2(19)); 
polyline([.3,xs),[.53,ys),'-1’,'sc’) 


c3=-32; 

s=[num2str(c3),' dB’); 
text(.2,.49,s,'sc’) 
{xs,ys]=dce2sc(PB(19),a3(19)); 


polyline, (.3,xs],[.50,ys],'-1’,'sc’) 
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C4x-29; 

s=(num2str(c4),' dB’); 
text(.2,.46,s,'sc’ 
(xs,ys}#dc2sc(PB(19),a4(19)); 
polyline([.3,xs),{.47,ys],'-r’,'sc’) 


c5=-26; 

s=([num2str(c5),' dB’); 
text(.2,.43,s,'sc’) 
{xs,ys]=de2sc(PB(19),a5(19)); 
polyline([.3,xs},[.44,ys],'-r’,'sc') 


c6=-23; 

s=[num2str(c6),’ dB’); 
text(.2,.40,s,'sc’) 
[xs,ys]=dc2sc(PB(19),a6(19)); 
polyline((.3,xs),{.41,ys],'-1’,'sc’) 


c7=-20; 

s=[num2str(c7),' dB’); 
text(.2,.37,s,'sc') 
[xs,ys}=de2sc(PB(19).27(19)); 


polyline([.3,xs},[.38,ys],'-r','sc’) 


meta t2as13xa 
pause 


axis([-60 -44 -21 -5)) 
plot(PB,p1,'-',PB,p2,'-',PB,p3,'-',PB,p4,'-',PB,p5,'-',PB,p6,'-’,PB,p7,°-’),.. 
semilogy, grid,. 

%otitle(‘case 1. 3,nu=. 1,delk=.3: Four channels-Asynchronous’),.. 
xlabel(‘Peak Power(dBW)' ),ylabel(‘Probability of Bit Error’) 


s=['Crosstalk='}; 
text(.15,.6,s,'sc’) 


cl=-20; 

s=[num2str(c1),’ dB‘); 
text(.2,.55,s,'sc') 
[xs,ys}=dce2sc(PB(15),p7(15)); 
polyline({.3,xs),[.56,ys),'-r’,'sc’) 


c2=-23; 

s=[num2str(c2),' dB’); 
text(.2,.52,s,'sc’) 
[xs,ys]=dc2sc(PB(15),p6(15)); 
polyline((.3,xs),[.53,ys),'-1’,'sc’) 


c3=-26; 

s=[num2str(c3),' dB’); 
text(.2,.49,s,'sc’) 
[xs,ys]=dce2sc(PB(15),p5(15)); 
polyline([.3,xs},[.50,ys],'-1’,'sc’) 








') 
(xs Nstedese(PB(15) pA(15)); 
polyline([.3,xs},[.47,ys],'-1,'sc’) 


c5=z-32; 

s=(num2srr(c5),' dB’); 
text(.2,.43,s,'sc’) 
{xs,ys]=dce2sc(PB(15),p3(15)); 
polyline([.3,xs),[.44,ys],'-1’,'sc’) 


c6=-35; 

s=(num2str(c6),’ dB’); 
text(.2,.40,s,'sc’) 
[xs,ys]=dce2sc(PB(15),p2(15)); 
polyline({.3,xs),[.41,ys],'-1,,'sc’) 


c7=0; 

s=[' 'num2str(c7)]; 
text(.2,.37,s,'sc’) 
[xs,ys]=dc2sc(PB(15),p1(15)); 
polyline([.3,xs],[-38,ys],'-r’,'sc’) 


axis; 
meta t2as13xb 











APPENDIX E - MATLAB MODEL FOR THREE ADJACENT 
CHANNELS WITH SYNCHRONOUS TRANSMISSION 


This model is for the case of three adjacent channels with (v,5,,82,5;) = 


(0.1,0.3,0.3,0.3) for synchronous transmission. The MATLAB version used was version 
rr as program itself is a translation of equations (1) through (13) into MATLAB code 
or k=3. 


% Synchronous Case 

% k=3 implying three channels 

% Case 1.3: nu=.1 and delta function subscript k (delk)=.3 
% th3ch13.m 

% 1 Aug 93 

% Thermal Noise=1E-22 AA2/Hz 


% Known constants 


R=.5; 
T=1/(S00E6); 
NO=1E-22; 
q=1.6E-19; 
Idk=10E-9; 


% Coupling constants for Ck*2 
ck1=0; 

ck2=3.16E-4; %-35 dB 
ck3=6.32E-4; %-32 dB 
ck4=1.26E-3; %-29 dB 
ck5=2.52E-3; %-26 dB 
ck5a=5.011E-3; %-23 dB 
ck6=.01; %-20dB 
ck7=1.995E-2; %-17 dB 
ck8=3.981E-2; %-14 dB 
ck9=7.943E-2; %-11 dB 
ck10=.1585; %-8 dB 

ck=[0 3.16E-4 6.32E-4 1.26E-3 2.52E-3 5.012E-3 .01 1.995E-2 3.981E-2 
7.943E-2 .1585]; 


% Pattern bit values 

%b1l=[.2182 .9759 .9759 .2182 .2182 .2182 .9759 .9759); 
%for b21,b11,b01 for patterns A,B,C,D,E,F,G,H 
%b2=[.9759 .2182 .9759 .2182 .2182 .9759 .2182 .9759); 
Gofor b22,b12,b02 for patterns A,B,C,D,E,F,G,H 
%b3=[.9759 .9759 .9759 .9759 .2182 .2182 .2182 .2182]); 
%for b23,b13,b03 

b1lA=.2182; 

b1B=.9759; 

b1C=.9759; 














b1D=.2182; 

blE=.2182; 

bl F=.2182; 

b1G=.9759; 

b]H=.9759; 

b2A=.9759; 

b2B=.2182; 

b2C=.9759; 

b2D=.2182; 

b2E=.2182; 

b2F=.9759; 

b2G=.2182; 

b2H=.9759; 

b3A=.9759; 

b3B=.9759; 

b3C=.9759, 

b3D=.9759; 

b3E=.2182; 

b3F=.2182; 

b3G=.2182; 

b3H=.2182; 
b4A=b1A42+b2A42+b3A‘2; 
b4B=b1B42+b2B42+b3B2; 
b4C=b1C42+b2C42+b3C42, 
b4D=b1D‘2+b2D%2+b3D“2; 
b4E=b1EA2+b2E%2+b3E“2; 
b4F=b1F42+b2F%2+b3F2, 
b4G=b1G42+b2G"2+b3G"2; 
b4H=b1H‘2+b2H“2+b3H“2; 


b00=.2182; % approximate zero 
b10=.9759; % approximate 1 


% Determine k constant for the variance of X 

nu=.1; 

delk=.3; 

ka=(2*pi*nu)/(4*pi42*(nu42+delk’2)); 
kb=1/(16*pi44*(nu42+delk*2)2); 

rn dE SN a A aca a 
2*pi*nu*cos(2*pi*delk)); 
kd=(2*pi*delk*exp(-2*pi*nu))*(2*pi*(delk*cos(2*pi*delk)+nu*sin(2*pi*delk))); 
ke=4*pi42*(nu*2-delk’2); 

k=ka-kb*(kc+kd+ke); 

n=]; 

=1; 

m=]; 


% -60dB<=Peak Power<=-44dB 
% Peak power is A“2/2 = Ps 
for ss=-60:.2:-44; 
ps(m)=10(ss/10); 
m=m-+1; 
end 
% The goal here is to solve the equation for the optimal threshold 
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% vs peak power 


ha=zeros(length(ps),length(ck)); 
PBE=zeros(length(ps),length(ck)); 


for j=1:length(ck); % coupling values loop 
n=]; 
for n=1:length(ps); % peak ppower values 
Ps=ps(n); 


% Constants 


qcnst=q/(Ps*T*F) 
Tnoise=NO/(RA2*Ps*2*T); 
vdark=(q*Idk)/(R42*Ps‘2*T); 


% Determine the sigma for each of the 0 bit patterns 


sOA=sqrt(4*ck(j)*b00*2*k*b4A+qcnst*(b00“2+ck(j)*b4A)+vnoise+vdark); 
celery, Shae PMgy abe reiDh gale sabes hor eb earrestfienslnt 
sqrt(4*ck(j)*b00*2*k*b4C+qcnst*(b00/2+ck(j)*b4C)+vnoise+vdark); 
Dean sqrt(4*ck(j)*b00*2*k*b4D+qcnst*(b0042+ck(j)*b4D)+vnoise+vdark); 
sOB=sart(4*ck 4)*b00*2*k*b4E+qnnst*(b0O*2+ck()*b4E)+vnoise+vdark), 
s0F=sqrt(4*ck(j)*b00*2*k*b4F+qcnst*(b00‘2+ck(j)*b4F)+vnoise+vdark); 
sOG=sqrt(4*ck(j)*b00*2*k*b4G+qcnst*(b00“2+ck(j)*b4G)+vnoise+vdark); 
sOH=sqrt(4*ck(j)*b00*2*k*B4H+qcenst*(b00“2+ck(j)*b4H)+vnoise+vdark); 


% Determine the variance for each of the 0 bit patterns 


SOA2=s0A*2; 


s0H2=s0H"2; 


% Determine the sigma for each of the 1 bit patterns 


s1A=sqrt(4*ck(j)*b1042*k*b4A+qcenst*(b1042+ck(j)*b4A)+vnoise+vdark); 
s1B=sqrt(4*ck(j)*b10°2*k*b4B+qcenst*(b10“2+ck(j)*b4B)+vnoise+vdark); 
$1C=sqrt(4*ck(j)*b10°2*k*b4C+qcenst*(b10*2+ck(j)*b4C)+vnoise+vdark); 
s1D=sqrt(4*ck(j)*b10°2*k*b4D+qenst*(b1042+ck(j)*b4D)+vnoise+vdark); 
$1E=sqrt(4*ck(j)*b10°2*k*b4E+qcnst*(b1042+ck(j)*b4E)+vnoise+vdark); 

s1F=sqrt(4*ck(j)*b10°2*k*b4F+qcenst*(b102+ck(j)*b4F)+vnoise+vdark); 

s1G=sqrt(4*ck(j)*b10°2*k*b4G+qenst*(b1042+ck(j)*b4G)+vnoise+vdark); 
s1H=sqrt(4*ck(j)*b1042*k*b4H+qenst*(b10°2+ck(j)*b4H)+vnoise+vdark); 


% Determine the variance for each of the 1 bit patterns 
s1A2=s1A42; 


s1B2=s1B42; 
s1C2=s1C2; 
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s1D2=s1D“2; 
sl1E2=slE‘2; 
s1F2=s1F42; 
$1G22s51G42; 
$31H2#s1H“2; 


% Determine the mean for each of the 0 bit patterns 
m0A=b00‘2+ck(j)*b4A; 
m0B=b0042+ck(j)*b4B; 
m0C=b00'2+ck(j)*b4C; 
m0D=b00"2+ck(j)*b4D; 
mOQE=b00/2+ck(j)*b4E; 
m0F=b00°2+ck(j)*b4F, 
m0G=b002+ck(j)*b4G; 
m0H=b0042+ck(j)*b4H; 





% Determine the mean for each of the 1 bit patterns 


m1A=b10*2+ck(j)*b4A; 
m1B=b1042+ck(j)*b4B; 
m1C=b10A2+ck(j)*b4C; 
m1D=b1042+ck(j)*b4D; 
mlE=b1042+ck(j)*b4E; 

m1F=b1042+ck(j)*b4F; 

m1G=b1042+ck(j)*b4G; 
m1H=b1042+ck(j)*b4H; 


%Solve the equation F(a)=0 to optimize a, the threshold. F(a)=0 comes 
% from (dP/da)=0 


a=0:.001:1; 
FA=(1/s0A)*exp(-((a-m0A).42)./(2*s0A2))-(1/s1A)*exp(-((m1A-a).42)/(2*s1A2)); 
FB=(1/sOB)*exp(-((a-m0B).‘2)./(2*sOB2))-(1/s1B)*exp(-((m1B-a).42)/(2*s1B2)); 
FC=(1/sOC)*exp(-((a-m0C).‘2)./(2*s0C2))-(1/s1C)*exp(-((m1C-a).42)/(2*s1C2)); 
FD=(1/sOD)*exp(-((a-m0D).42)./(2*sOD2))-(1/s1D)*exp(-((m1D-a).42)/(2*s1D2)); 
FE=(1/s0E)*exp(-((a-m0E).*2)./(2*sOE2))-(1/s1E)*exp(-((m1E-a).42)/(2*s1E2)); 

=(1/sOF)*exp(-((a-m0F).42)./(2*sOF2))-(1/s1F)*exp(-((m1F-a).“2)/(2*s1F2)); 
FG=(1/s0G)*exp(-((a-m0G).42)./(2*s0G2))-(1/s1G)*exp(-((m1G-a).42)/(2*s1G2)); 
FH=(1/sOH)*exp(-((a-m0H).42)./(2*s0H2))-(1/s1H)*exp(-((m1H-a).42)/(2*s1H2)); 
F=FA+FB+FC+FD+FE+FF+FG+FH; 


[Y,1}=min(abs(F)); 
alpha(n,j)=a(); 
aa=alpha(n,j); 
% Determine the probability of bit error for each optimal threshold 
% over all bit patterns A,B,C,D 
pOA(n)=erf((alpha(n,j)-m0A)/sOA. inf); 
pOB(n)=erf((alpha(n,j)-m0B)/sOB, inf); 
pOC(n)=erf((alpha(n,j)-m0C)/sOC, inf); 
pOD(n)=erf((alpha(n,j)-m0D)/sOD, inf); 
pOE(n)=erf((alpha(n,j)-m0E)/sOE, inf); 
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pOF(n)=erf((alpha(n,j)-mOF)/sOF inf); 
pOG(n)=erf((alpha(n,j)-m0G)/sOG, inf); 
pOH(n)=erf((alpha(n,j)-m0H)/sOH, inf); 


plA(n)=erf((m1A-alpha(n,j))/s1A,inf); 
p1B(n)=erf((m1B-alpha(n,j)/s1B,inf); 
p1C(n)=erf((m1C-alpha(n,j))/s1C, inf); 
p1D(n)=erf((m1D-alpha(n,j)/s1D,inf); 
plE(n)=erf((m1E-alpha(n,j))/s1E,inf); 

p1F(n)=erf((m1F-alpha(n,j)Vs1F,inf); 

p1G(n)=erf((m1G-alpha(n,j))/s1G, inf); 
p1H(n)=erf((m1H-alpha(n,}))/s1H, inf); 


PBE(n,j)=(pOA(n)+p1A(n)+p0B(n)+p1B(n)+p0C(n)+p1C(n)+p0D(n)+p1D(n) 
+p0E(n)+p1E(n)+ EAU RAM LR rs Hy pHa ‘ 


end % Ps loop 
end %Cxk loop 


al=alpha(:,1); 
a2=alpha(:,2); 
a3=alpha(:,3); 
a4=alpha(:,4); 
a5=alpha(:,5); 
a6=alpha(:,6); 
a7=alpha(:,7); 
a8=alpha(:,8); 
a9=alpha(:,9); 
alQ=alpha(:,10); 
all=alpha(:,11); 


p1=(PBE(:,1)); 
p2=(PBE(:,2)); 
p3=(PBE(:,3)); 
p4=(PBE(:,4)); 
p5=(PBE(C,5)); 
p6=(PBE(:,6)); 
p7=(PBE(:,7)); 
p8=(PBE(:,8)); 
p9=(PBE(:,9)); 
p10=PBE(:,10); 
pl1=PBE(:,11); 


PB=10*log10(ps); 
plot(PB,al,'-',PB,a2,'-',PB,a3,'-',PB,a4,'-',PB,a5,'-',PB,a6,'-',PB,a7,'~’),.. 


grid... 
Stitle(case 1.3,nu=.1,delk=.3: Three channels’)... 
xlabel(‘Peak Power(dBW)'),ylabel(Optimal Threshold’) 


=('Crosstalk='}; 
text(.15,.6,s,'sc’) 
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cl=0, 

s=[" ‘,num2str(cl)); 
text(.2,.55,S,'sc’) 
[xs,ys}odedsc(PB(48),a1(48)); 
polyline([.3,xs},[.56,ys),'-r’,'sc’) 


c2=-35; 

s=(num2str(c2),’ dB’); 
text(.2,.52,s,'sc') 
[xs,ys]=dc2sc(PB(48),a2(48)); 
polyline([.3,xs],[.53,ys],'-r’,'sc’) 


c3=-32; 

s=[{num2str(c3),' dB’); 
text(.2,.49,s,'sc’) 
[xs,ys]=dc2sc(PB(48),a3(48)); 
polyline([.3,xs],[.5,ys],’-r’,'sc’) 


c4=-79; 

s=[num2str(c4),' dB’); 
text(.2,.46,s,'sc’) 

[xs, ys]=dc2sc(PB(48),a4(48)); 
polyline([.3,xs},[.47,ys],'-r’,'sc’) 


c5=-26 

sainummbeieo): dB’); 
text(.2,.43,s,'sc’) 
[xs,ys]=dc2sc(PB(48),a5(48)); 
polyline([.3,xs},[.44,ys],'-r’,'sc’) 


c6=-23; 

s=(num2str(c6), dB'); 
text(.2,.4,s,'sc’) 
[xs,ys]ade2se(PB(48), a6(48)); 


polyline((.3,xs],[.41,ys],'-1','sc’) 


c7=-20; 

s=[num2str(c7),' dB’); 
text(.2,.37,s,'sc') 
[xs,ys}=dc2sc(PB(48),a7(48)); 
acahrpe aa .3,xS},[.38,ys],'-r','sc’) 


pause 

axis([-60 -44 -21 -5)) 
plot(PB,p1,'-',PB,p2,'-',PB,p3,'-',PB,p4,'-',PB,p5,'-',PB,p6,'-',PB,p7,'-'),.. 
semilogy,grid, 

%title(‘case 1 O,nu=. 1,delk=0: Three channe)s’)... 
xlabel(Peak Power(dBW)’),ylabel( Probability of Bit Error’) 


s=['Crosstalk='J; 
text(.15,.6,s,'sc’) 
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cl=-20; 

s=[num2str(cl),’ dB’); 
text(.2,.55,s,'sc’) 
{xs,ys}=dc2sc(PB(38),p7(38)); 
polyline({.3,xs},[.56,ys],'-1','sc’) 


c2=-23; 

s=[num2str(c2),' dB’]; 
text(.2,.52,s,'sc’) 
[xs,ys}=dc2sc(PB(38),p6(38)); 
polyline((.3,xs],{.53,ys],'-r',’sc’) 


c3=-26; 

s=[num2st(c3),’ dB’); 
text(.2,.49,s,'sc') 
(xs,ys]=dc2sc(PB(38),p5(38)); 
polyline([.3,xs],{.5,ys],'-1','sc’) 


c4=-29; 

s=[num2str(c4),’ dB’); 
text(.2,.46,s,'sc’) 
[xs,ys]=dc2sc(PB(38),p4(38)); 
polyline((.3,xs],[.47,ys],'-1','sc’) 


c5=-32; 

s=[num2str(c5),' dB’); 
text(.2,.43,s,'sc’) 
[xs,ys]=dc2sc(PB(38),p3(38)); 
polyline({.3,xs],[.44,ys],'-r','sc’) 


c6=-35; 

s=[num2str(c6),’ dB’); 
text(.2,.40,s,'sc’) 
[xs,ys]=dc2sc(PB(38),p2(38)); 
polyline([.3,xs],[.41,ys],'-1','sc’) 


c7=0; 

s=[' ‘num2str(c7)]; 
text(.2,.37,s,'sc') 
[xs,ys]=dc2sc(PB(38),p1(38)); 
polyline({.3,xs],[.38,ys],'-1','sc’) 
axis; 


meta t3ch13xb 
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APPENDIX F - MATLAB MODEL FOR FOUR ADJACENT 
CHANNELS FOR SYNCHRONOUS TRANSMISSION 


This model is for the case of four adjacent channels with (v,5,,52,53,8,) = 


(0.1,0.3,0.3,0.3,0.3) for synchronous transmission. The MATLAB version used was 
version 3.5k. The program itself is a translation of equations (1) through (13) into 
MATLAB code for k=4. 





% Synchronous Case 

% k=4 implying four channels 

% Case 1.3: nu=.1 and delta function subscript k (delk)=.3 
% th4ch13x.m 

% 22 Aug 93 

% Thermal Noise: 1E-22 A42/Hz 


% Known constants 


R=.5; 
T=1/(500E6); 
NO=1E-22; 
q=1.6E-19; 
Idk=10E-9; 


% Coupling constants for Ck42 


ck2=3.16E-4; %-35 dB 
ck3=6.32E-4; %-32 dB 
ck4=1.26E-3; %-29 dB 
ck5=2.52E-3; %-26 dB 
ckSa=5.011E-3; %-23 dB 
ck6=.01; %-20dB 
ck7=1.995E-2; %-17 dB 
ck8=3.981E-2; %-14 dB 
ck9=7.943E-2; %-11 dB 
ck10=.1585; %-8 dB 
ck=[0 3.16E-4 6.32E-4 1.26E-3 2.52E-3 5.011E-3 .01 1.995E-2 3.981E-2 
7.943E-2 .1585); 


% Pattern bit values 

%b1=[.2182 .9759 .9759 .2182 .2182 .2182 .9759 9759 .2182 .2182 .2182 
.9759 .9759 .9759 .2182 .9759]; 

%for b31,b21,b11,b01 for patterns A,B,C,D,E,F,G,H,I,J,K,L,M,N,O,P 

Gob2=[.9759 2182 9759 .2182 .2182 .9759 .2182 .9759 .2182 .2182 .9759 
.2182 .9759 .2182 .9759 .9759]); 

%for b32,b22,b12,b02 for patterns A,B,C,D,E,F,G,H,1,J,K,L,M,N,O,P 

%b3=[.9759 .9759 .9759 .9759 .2182 .2182 .2182 .2182 .2182 .9759 .2182 
.2182 .2182 .9759 .9759 .9759]; 

for b33,523,b13,b03 
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%b4=(.9759 9759 .9759 .9759 .9759 9759 .9759 9759 .2182 .2182 .2182 
.2182 .2182 .2182 .2182 .2182}; 
%for b34,b24,b14,b04 
b1A=.2182; 
b1B=.9759; 
b]1C=.9759; 
b1D=.2182; 
b1lE=.2182; 
b1F=.2182; 
b1G=.9759; 
b1H=.9759; 
blI=.2182; 
b1J=.2182; 
b1K=.2182; 
b1L=.9759; 
b1M=.9759; 
b1IN=.9759; 
b10=.2182; 
b1P=.9759; 
b2A=.9759; 
b2B=.2182; 
b2C=.9759; 
b2D=.2182; 
b2E=.2182; 
b2F=.9759; 
b2G=.2182; 
b2H=.9759; 
b2]=.2182; 
b2J=.2182; 
b2K=.9759; 
b2L=.2182; 
b2M=.9759; 
b2N=.2182; 
b20=.9759; 
b2P=.9759; 
b3A=.9759; 
b3B=.9759; 
b3C=.9759; 
b3D=.9759; 
b3E=.2182; 
b3F=.2182; 
b3G=.2182; 
b3H=.2182; 
b3I=.2182; 
b3J=.9759; 
b3K=.2182; 
b3L=.2182; 
b3M=.2182; 
b3N=.9759; 
b30=.9759; 
b3P=.9759; 
b4A=.9759; 
b4B=.9759; 
b4C=.9759; 
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b4D=.9759, 
b4E=.9759, 
b4F=.9759; 
b4G=.9759; 
b4H=. 9759; 


b4P=. 2182: 


b5A=b1A42+b2A42+b3A42+b4A‘2; 
b5B=b1B42+b2B42+b3B42+b4B42; 
bSC=b1C42+b2C42+b3C42+b4C42; 
b5D=b1D42+b2D42+b3D‘2+b4D“2; 
bSE=b1E42+b2E%2+b3E%2+b4E‘2; 
bSF=b1F42+b2F42+b3F42+b4F42; 
b5G=b1G42+b2G42+b3G42+b4G42; 
bSH=b1H42+b2H‘2+b3H‘2+b4H42; 
b51=b1H42+b2H42+b3H42+b4142; 
b5J=b1JA2+b2J42+b3IA2+b4)2; 
b5K=b1K42+b2K%2+b3K‘2+b4K‘2; 
b5L=b1L42+b2L42+b3L“2+b4L‘2; 
bSM=b1M‘42+b2M‘42+b3M‘2+b4M‘2; 
bSN=b1N42+b2N‘42+b3N‘42+b4N‘2; 
b5O=b1042+b20%2+b30%2+b40°2; 
b5P=b1P42+b2P42+b3P42+b4P"2; 


b00=.2182; % approximate zero 
b10=.9759; % approximate 1 


% Determine k constant for the variance of X 

nu=.1; 

delk=.3; 

ka=(2*pi*nu)/(4*pi42*(nu42+delk“2)); 

kb=1/(16*pi44*(nu42+delk*2)2); 
kce=(2*pi*nu*exp(-2*pi*nu))*(2*pi*delk*sin(2*p i*delk)-2*pi*nu*cos(2*pi*delk)); 
kd=(2*pi*delk*exp(-2*pi*nu))*(2*pi*(delk*cost2*p i*delk)+nu*sin(2*pi*delk))); 
ke=4*pi42*(nu42-delk’2); 

k=ka-kb*(kc+kd+ke); 

n=]; 

1; 

m=1; 


% -60dB<=Peak Power<=-44dB 
% Peak power is A‘2/2 = Ps 
for ss=-60:.5:-44; 
ps(m,- 10(ss/10); 
m=m-+1; 
end 
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% The goal here is to solve the equation for the optimal threshold 
% vs peak power 


alpha=zeros(length(ps),length(ck)); 
PBE=zeros(length(ps),length(ck)); 


for I slength(ck); % coupling values loop 
net: 


for n=! slength(ps); % peak ppower values 
Ps=ps(n); 
% Constants 


qcnst=q/(Ps*T*R); 
vnoise=N0/(R“2*Ps*2*T); 
vdark=(q*Idk)/(R“2*Ps42*T); 


% Determine the sigma for each of the 0 bit patterns 


sOA=sqrt(4*ck(j)*b00°2*k*b5A+qcenst*(b00‘2+ck(j)*b5A)+vnoise+vdark); 
sOB=sqrt(4*ck(j)*b0042*k*b5B+qcnst*(b00“2+ck(j)*b5B)+vnoise+vdark); 
s0C=sqrt(4*ck(j)*b00*2*k*b5C+qcnst*(b00.2+ck(j)*b5C)+vnoise+vdark); 
sOD=sqrt(4*ck(j)*b0042*k*b5D+qcnst*(b00"2+ck(j)*b5D)+vnoise+vdark); 
sOE=sqrt(4*ck(j)*b00A2*k*bS5E+qcnst*(b00“2+ck(j)*bSE)+vnoise+vdark); 
sOF=sqrt(4*ck(j)*b00A2*k*b5SF+qcnst*(b00“2+ck(j)*b5F)+vnoise+vdark); 
s0G=sqrt(4*ck(j)*b00*2*k*b5G+qenst* (b00“2+ck(j)*b5G)+vnoise+vdark); 
sOH=sqrt(4*ck(j)*b0042*k*bSH+qcenst*(b00%2+ck(j)*b5H)+vnoise+vdark); 
sOI=sqrt(4*ck(j)*b00*2*k*bSI+qcnst*(b0042+ck(j)*b51)+vnoise+vdark); 
sOJ=sqrt(4*ck(j)*b00*2*k*b5J+qenst*(b0042+ck(j)*b5J)+vnoise+vdark); 
sOK=sqrt(4*ck(j)*b00*2*k*b5K+qcenst*(b0042+ck(j)*bSK)+vnoise+vdark); 
sOL=sqrt(4*ck(j)*b00‘2*k*b5L+qcnst*(b0042+ck(j)*b5L)+vnoise+vdark); 
sOM=sqrt(4*ck(j)*b00*2*k*b5M+qcens. *(60042+ck(j)*b5M)+vnoise+vdark); 
sON=sqrt(4*ck(j)*b00“2*k*b5N+qcnst*(b0042+ck(j)*b5N)+vnoise+vdark); 
s0O=sqrt(4*ck(j)*b00*2*k*b50+qenst*(b00*2+ck(j)*b50)+vnoise+vdark); 
sOP=sqrt(4*ck(j)*b00*2*k*b5P+qcenst*(b002+ck(j)*b5P)+vnoise+vdark); 


% Determine the variance for each of the 0 bit patterns 


sOA2=s0A4‘2; 
s0B2=s0B“2; 
s0D2=s0D“2; 
sOE2=sO0E“2; 
sOF2=sOF‘2; 
s0G2=s0G“2; 
sOH2=s0H“2; 
sQI2=sO0I42; 
s0J2=s0J‘2; 
sOK2=s0K‘2; 
sOL2=s0L‘2; 
s0M2=s0M“2; 
sON2=sON42; 
s002=s00“2; 
sOP2=s0P‘2; 
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% Determine the sigma for each of the 1 bit patterns 


$1 A=sqrt(4*ck(j)*b1042*k*b5A+qcnst*(b1042+ck(j)*b5A)+vnoise+vdark); 
$1 B=sqrt(4*ck(})*b1042*k*b5B+qcenst*(b10°2+ck(j)*b5B)+vnoise+vdark); 
s1Cxsqrt(4*ck(j)*b10°2*k*b5C+qcnst*(b1042+ck(j)*b5C)+vnoise+vdark); 
$1D=sqrt(4*ck(Q)*b10°2*k*b5D+qcnst*(b1042+ck(j)*b5D)+vnoise+vdark); 
$1 E=sqrt(4*ck(j)*b10°2*k*bSE+qcnst*(b]0%2+ck(})*bSE)+vnoise+vdark); 
siFe=sqrt(4*ck(j)*b10°2*k*bSF+qcnst*(b10%2+ck(j)*b5F)+vnoise+vdark); 
$1G=sqrt(4*ck(j)*b10°2*k*b5G-+qcnst*(b10\2+ck(j)*b5G)+vnoise+vdark); 
s1H=sqrt(4*ck(j)*b10°2*k*b5H+qcenst*(b1042+ck(j)*b5H)+vnoise+vdark); 
slI=sqrt(4*ck(j)*b10°2*k*b51+qcnst*(b10°2+ck(j)*b51)+vnoise+vdark); 
slJ=sqrt(4*ck(j)*b10*2*k*b5J+qcenst*(b1042+ck(j)*b5J)+vnoise+vdark); 
slK=sqrt(4*ck(j)*b1042*k*b5K-+qcenst*(b1042+ck(j)*b5K)+vnoise+vdark); 
s1]L=sqrt(4*ck(j)*b1042*k*bSL+qcnst*(b10%2+ck(j)*bSL)+vnoise+vdark); 
s1M=sqrt(4*ck(j)*b10°2*k*b5M-+qcnst*(b10%2+ck(j)*b5M)+vnoise+vdark); 
S1N=sqrt(4*ck(j)*b10°2*k*b5N-+qcnst*(b1042+ck(j)*b5N)+vnoise+vdark); 
s1O=sqrt(4*ck(j)*b1042*k*b50+qcnst*(b1042+ck(j)*b50O)+vnoise+vdark); 
$1P=sqrt(4*ck(j)*b10°2*k*b5P+qcnst*(b1042+ck(j)*b5P)+vnoise+vdark); 


% Determine the variance for each of the 1 bit patterns 


sl1A2=s1A42; 
$1B2=s1B42; 
s1C2=s1CA2; 
s1D2=s1D/2; 
slE2=slE‘2; 
s1F2=s1F42; 
$1G2=s1G“2; 
S1H2=s1H‘2; 
slI2=s1I42; 
slJ2=s1J42; 
s1K2=s1K‘2; 
s1L2=s1L‘2; 
$1M2=s1M‘2; 
sIN2=s1N42; 
s102=s10/2; 
s1P2=s1P42; 


% Determine the mean for each of the 0 bit patterns 


mOA=b0012+ck(j)*b5A; 
m0B=b002+ck(j)*b5B; 
m0C=b0012+ck(j)*b5C; 
m0D=b00\2+ck(j)*bSD; 
mOE=b0042+ck(j)*b5E; 
m0F=b00*2+ck(j)*b5F; 
mOG=b00*2+ck(j)*b5G; 


m0J=b0042+ck(j)*b5J; 
m0K=b002+ck(j)*b5K; 
m0L=b00°2+ck(j)*b5L; 
m0M=b00\2+ck(j)*b5M:; 
mON=b0012+ck(j)*bSN; 
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m00=b00*2+ck(j)*b50; 
mOP=b0042+ck(j)*bSP; 


% Determine the mean for each of the 1 bit patterns 


mlA=b10*2+ck(j)*bSA; 
m1B=b10%2+ck(j)*b5B; 
m1C=b1042+ck(j)*b5C; 
m1D=b10°2+ck(j)*b5D; 
mlE=b102+ck(j)*bSE; 
m1F#b1042+ck(j)*b5F; 
m1G=b10*2+ck(j)*b5G; 
m1H=b1042+ck(j)*b5H; 
mlI=b10*2+ck(j)*b5I; 
m1lJ=b1042+ck(j)*b5J; 
m1K=b1042+ck(j)*b5K; 
m1L=b10°2+ck(j)*bSL; 
m1M=b102+ck(j)*b5M; 
m1N=b1042+ck(j)*v5N; 
m10=b]042+ck(j)*b5O; 
m1P=b10°2+ck(j)*b5P; 


%Solve the equation F(a)=0 to optimize a, the threshold. F(a}=0 comes 
% from (dP/da)=0 


a=0:.0009:1; 
FA=(1/s0A)*exp(-((a-m0A).42)./(2*s0A2))-(1/s1A)*exp(-((m1A-a).42)/(2*s1A2)); 
FB=(1/sOB)*exp(-((a-m0B).42)./(2*s0B2))-(1/s1B)*exp(-((m1B-a).42)/(2*s1B2)); 
FC=(1/s0C)*exp(-((a-m0C).42)./(2*s0C2))-(1/s1C)*exp(-((m1C-a).42)/(2*s1C2)); 
FD=(1/sOD)*exp(-((a-m0D).‘2)./(2*s0D2))-(1/s1D)*exp(-((m1D-a).42)/(2*s1D2)); 
FE=(1/sOE)*exp(-((a-m0E).*2)./(2*sOE2))-(1/s1E)*exp(-((m1E-a).42)/(2*s1E2)); 

=(1/sOF)*exp(-((a-mOF).42)./(2*sOF2))-(1/s1F)*exp(-((m1F-a).42)/(2*s1F2)); 
FG=(1/s0G)*exp(-((a-m0G).“2)./(2*s0G2))-(1/s1G)*exp(-((m1G- a).42)/(2*s1G2)); 
FH=(1/s0H)*exp(-((a-m0H).‘2)./(2*sOH2))-(1/s1H)*exp(-((m1H-a).42)/(2*s1H2)); 
Fi=(1/sOI)*exp(-((a-mOI).42)./(2*s012))-(1/s11)*exp(-((m1I-a).42)/(2*s112)); 
FJ=(1/s0J)*exp(-((a-m0J).42)./(2*s0J2))-(1/s1J)*exp(-((m1J-a).42)/(2*s1J2)); 
FK=(1/sOK)*exp(-((a-mOK).42)./(2*sOK2))-(1/s1K)*exp(-((m1K-a).42)/(2*s1K2)); 
FL=(1/sOL)*exp(-((a-m0L).42)./(2*sOL2))-(1/s1L)*exp(-((m1L-a).42)/(2*s1L2)); 
FM=(1/sOM)*exp(-((a-m0M).42)./(2*s0M2))-(1/s1M)*exp(-((m1M- 
a).2)/(2*s1M2)); 
FN=(1/sON)*exp(-((a-mON).42)./(2*sON2))-(1/s1N)*exp(-((m1N-a).42)/(2*s1N2)); 
FO=(1/s0O)*exp(-((a-m00).42)./(2*s002))-(1/s10)*exp(-((m10-a).42)/(2*s102)); 
FP=(1/sOP)*exp(-((a-mOP).42)./(2*sOP2))-(1/s1P)*exp(-((m1P-a).42)/(2*s1P2)); 
F=FA+FB+FC+FD+FE+FF+FG+FH+FI+FJ+FK+FL+FM+FN+FO+FP; 


[Y,I]}=min(abs(F)); 
alpha(n,j)=a(); 
aa=alpha(n,j); 


% Determine the probability of bit error for each optimal threshold 
% over all bit pattems A,B,C,D 
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pOA (n)=erf((alpha(n,j)-m0A )/sOA inf); 
pOB(n)=erf((alpha(nj)-m0B)/sOB, inf); 
pOC(n)=erf((alpha(n,})-m0C)/sOC inf); 
pOD(n)=erf((alpha(n,j)-m0D)/sOD, inf); 
pOE(n)=erf((alpha(n,j)-m0E)/sOE, inf); 
pOF(n)=erf((alpha(n,j)-m0F)/sOF, inf); 
pOG(n)=erf((alpha(n,j)-m0G)/s0G, inf); 
Ceti aps Ui slr 
pOl(n)=erf((alpha(n,j)-m0I)/sOI,inf); 
p0J(n)=erf(( re ay per dane 
POR en pha(n,j)-m0K)/sOK, inf); 
pOL(n)=erf((al 2 arene ay ip ee 
POM (pearl pha(n,j)-m0M)/sOM, inf); 
PON pect iatonaGnyONSON inf); 
p0O(n)=erf((alpha(n,j)-m00)/s0O, inf); 
pOP(n)=erf((alpha(n,j)-mOP)/sOP, inf); 


p1A(n)=erf((m1A-alpha(n,j))/s1A, inf); 
p1B(n)=erf((m1B-alpha(n j))/s1B inf); 
p1C(n)=erf((m1C-alpha(n,j))/s1C,inf); 
p1D(n)=erf((m1D-alpha(n,j))/s1D, inf); 
p1E(n)=erf((m1E-alpha(n,j))/s1E,inf); 
p1F(n)=erf((m1F-alphanj))/s1F,irf); 
p1G(n)=erf((m1G-alpha(n,j))/s1G,inf); 
p1H(n)=erf((m1H-alpha(n,j))/s1H.inf); 
p11(n)=erf((m1I-alpha(n,j))/s11,inf); 
p1J(n)=erf((m1J-alpha(n,j))/s1J,inf); 
p1K(n)=erf((m1K-alpha(n,j))/s1K,inf); 
p1L(n)=erf((m1L-alpha(n,j))/s1L,inf); 
p1M(n)=erf((m1M-alpha(nj))/s1M.inf); 
PIN(n)=erf((m1N-alpha(n,))/sIN sinf); 
p10(n)=erf((m10-alpha(n,j))/s10, inf); 
p1P(n)=erf((m1P-alpha(n,j))/s1P,inf); 


pA(n)=p0A(n)+p1A(n); 
pB(n)=p0B(n)+p1B(n); 
pC(n)=p0C(n)+p1C(n); 
pD(n)=p0D(n)+p1D(n); 
pE(n)=p0E(n)+p1E(n); 

pF(n)=p0F(n)+p1F(n); 
pG(n)=p0G(n)+p1G(n); 
pH(n)=p0H(n)+p1H(n); 
pl(n)=p0l(n)+p11(n); 

pJ(n)=p0J(n)+p1J(n); 

DAO aa 
pL(n)=p0L(n)+p1L{n); 
PM(n)=p0M(n)+p1M(n); 
pO(n)=p0O(n)+p10(n): 
pP(n)=p0P(n)+p1P(n); 








PEC 0) EE PC Bee a a Pome erpne (n) 
+pK(n)+pL(n)+pM(n)+pN(n)+pO(n)+pP(n))/64; 











end % Ps loop 
end %Ck loop 


al=alpha(: Ds 


a6=alpha(: °6); 
a7=alpha(: 7); 


a10=aiphac: 10); 
all=alpha(:,11); 


pl=(PBE(:,1)); 
p2=(PBE(:,2)); 
p3=(PBE(:,3)); 
p4=(PBE(:,4)); 
p5=(PBE(:,5)); 
p6=(PBE(:,6)); 
p7=(PBE(:,7)); 
p8=(PBE(: 8)); 
p9=(PBE( 9)); 
P10=PBEC 10); 
pll=PBE(:,11); 


PB=10*log10(ps); 
plot(PB,al .'-',PB,a2,'-',PB,a3,'-',PB,a4,'-',PB,a5,'-’,PB,a6,'-",PB,a7,'-'),.. 


grid... 
Sotitle(‘case 1.3,nu=.1,delk=.3: Four channels’)... 
xlabel(‘Peak Power(dBW)'),ylabel(‘Optimal Threshold’) 


s=[‘Crosstalk=']; 
text(.15,.6,s,'sc’) 


cl=0; 

s=[' ‘,ynum2str(c1)); 
text(.2,.55,s,'sc') 

[xs, ys]=de2sc(PB(19),a1 (19)); 


nolyline([. 3,xs],[.56,ys],'-r','sc’) 


c2=-35; 

s=[num2str(c2),’ dB'); 
text(.2,.52,s,'sc’) 
[xs,ys]=de2sc(PB(19),a2(19)); 
polyline([.3,xs],[.53,ys],'-r','sc’) 


c3=-32; 


s=[num2str(c3),' dB’); 
text(.2,.49,s,'sc’) 
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(xs,ys}=dc2sc(PB(19),a3(19)); 
polyline([.3,xs),{.50,ys],'-r’,'sc’) 


c4=-29; 

s=(num?str(c4),’ dB’); 
text(.2,.4€,s,'sc’) 
[xs,ys]=acZsc(PB(19),a4(19)); 


polyline({.3,xs],[.47,ys],'-1','sc’) 


c5=-26; 

s=[num2str(c5),' dB’); | 
text(.2,.43,s,'sc’) 

[xs,ys]=dc2sc(PB(19),a5(19)); 

polyline([.3,xs],[.44,ys],'-r’,'sc’) 


c6=-23; 

s=[num2str(c6),’ dB’); 
text(.2,.40,s,'sc’) 
{xs,ys]=dc2sc(PB(19),a6(19)); 


polyline({.3,xs},[.41,ys),'-1','sc’) 


c7=-20; 

s=[num2str(c7),' dB’); 
text(.2,.37,s,'sc') 
[xs,ys]=dc2sc(PB(19),a7(19)); 
polyline([.3,xs],[.38,ys],'-r','sc’) 


meta t4ch13xa 
pause 


axis([-60 -44 -21 -5]) 

plot(PB,p1,'-',PB,p2,'-',PB,p3,'-',PB,p4,'-',PB »p5,'-',PB,p6,'-',PB,p7,'-'),.. 
%PB,p8,'-',PB,p9,'-',PB,p10,'-',PB,p11,'-') 

semilogy,grid,.. 

%title(‘case 1.3,nu=.1,delk=.3: Four channels’)... 

xlabe](Peak Power(dBW)'),ylabel(‘Probability of Bit Error’) 


s=('Crosstalk='J; 
text(.15,.6,s,'sc’) 


cl=-20; 

s=(num2str(cl),’ dB’); 
text(.2,.55,s,'sc') 
[xs,ys]=dc2sc(PB(15),p7(15)); 
polyline({.3,xs],{.56,ys],'-r','sc’) 


c2=-23; 

s=[num2str(c2),' dB’); 
text(.2,.52,s,'sc’) 
[xs,ys}=dce2sc(PB(15),p6(15)); 
polyline((.3,xs],[.53,ys],'-1','sc’) 


c3=-26; 
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s=(num2str(c3),’ dB’); 
text(.2,.49,s,'sc 
[xs,ys}=de2sc(PB(15),p5(15)); 
polyline((.3,xs},[.50,ys],'-r’,'sc’) 


c4=-29; 

s=(num2str(c4),' dB’); 
text(.2,.46,s,'sc’) 
[xs,ys]=dc2sc(PB(15),p4(15)); 
polyline([.3,xs],[.47,ys],'-r',’sc’) 


c5=-32; 

s=[num2str(c5),’ dB’); 
text(.2,.43,s,'sc') 
[xs.ys]=dodse(PB(15),p3(15)); 
polyline([.3,xs],[.44,ys],'-1','sc’) 


c6=-35; 

s=[num2str(c6),’ dB’); 
text(.2,.40,s,'sc’) 
({xs,ys]=dc2sc(PB(15),p2(15)); 
polyline({.3,xs],[.41,ys],'-1','sc’) 


c7=0, 

s=[' ‘ num2str(c7)}; 
text(.2,.37,s,'sc’) 

[xs, yslede3sc(PB(15) »p1(15)); 


polyline({. 3,xs],[.38,ys],'-r','sc’) 
aie t4ch13xb 
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